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FOREWORD 


In  this  paper  computational  procedures  that  solve  three  general 
mathematical  problems  concerning  optimal  flows  through  networks  are 
described.  The  presentation  is  directed  toward  readers  interested  in 
complete  mathematical  validation  of  these  algorithms  or  in  the 
mathematical  details  of  their  implementation. 

The  requirement  for  incorporation  of  such  procedures  arose  in  the 
context  of  RAC  Study  24,1,  and  the  work  reported  here  represents  one 
aspect  of  that  study.  The  objective  of  the  study  is  to  design  a  computer- 
assisted  planning  procedure  for  determining  US  Army  midrange  require¬ 
ments  for  technical-service  troops  and  class  IV  materiel  for  specific 
campaigns.  One  of  the  tasks  involved  in  the  planning  process  is  the 
detailed  routing  of  supplies  through  highway  and  railroad  networks.  The 
material  presented  here  provides  a  mathematical  basis  for  efficiently 
selecting  desirable  routings  through  such  networks. 

The  first  of  the  three  procedures  determines  minimum-cost  flows 
between  specified  network  nodes  when  a  maximum  flow  rate  or  "  capacity  " 
is  given  for  each  link,  and  when  the  cost  of  flow  in  a  link  is  proportional 
to  the  magnitude  of  flow.  The  second  determines  minimum-cost  flows 

iii 


when  minimum  as  well  as  maximum  flow  rates  are  stipulated  for  all 
links.  The  third  determines  whether  or  not  a  given  flow  pattern  minimizes 
cost.  In  the  event  it  does  not,  the  procedure  then  transforms  the  given 
flow  pattern  into  one  that  is  optimal. 

A  computer  program  employing  the  first  two  of  these  procedures  has 
been  written  for  the  IBM  7090  computer  and  constitutes  one  operational 
package  in  the  new  planning  system.  Readers  interested  in  the  appli¬ 
cations  of  these  results  in  the  specific  context  of  Study  24.1  are  referred 
to  ORO-T-393,  "  Computer-Assisted  Strategic  Logistic  Planning;  Trans¬ 
portation  Phase";  ORO-SP-160,  "  Programming  Computer-Assisted 
Transportation  Planning";  and  ORO-TP-45,  "  Computer-Assisted  Prepa¬ 
ration  of  Transportation  Annex,  Planner  Manual.  "  Further  discussion, 
including  documentation  and  use  of  the  computer  program,  will  appear  in 
the  final  report  of  Study  24.1. 

The  material  presented  here  constitutes  an  adaptation  and  unified 
presentation  of  various  known  results  in  network  flow  theory,  in  a  form 
useful  to  support  the  above  study. 
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THREE  GENERAL  NETWORK  FLOW 
PROBLEMS  AND  THEIR  SOLUTIONS 


ABSTRACT 


Computational  procedures  for  solving  three  general  network  flow 
problems  are  presented,  together  with  proofs  establishing  their  validity. 
Two  of  the  problems  are  concerned  with  the  determination  of  feasible 
flows  (i.e. ,  flows  that  lie  between  prescribed  bounds  in  every  arc  of  the 
network)  whose  costs  are  minimum,  where  arc  costs  are  proportional 
to  the  magnitude  of  arc  flow.  The  third  problem  involves  the  determi¬ 
nation  of  whether  or  not  a  given  feasible  flow  has  minimum  cost.  The 
utility  of  the  three  procedures  is  illustrated  in  the  context  of  a  general 
transportation  application.  The  material  presented  constitutes  an  adap¬ 
tation  and  unified  discussion  of  various  known  concepts  and  results  in 
network  flow  theory. 
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INTRODUCTION 


This  paper  extends  the  methods  described  in  ORO-TP-15^  to  a  wider 
class  of  network  flow  problems.  In  that  publication,  networks  that  had  a 
specified  "  capacity"  as  well  as  a  "  unit  cost"  associated  with  each  direc¬ 
tion  of  every  edge  or  "link"  were  considered.  For  convenience,  parallel 
links  were  excluded  in  order  that  a  link  might  be  uniquely  identified  by  its 
two  vertices.  The  problem  under  consideration  concerned  the  determi¬ 
nation  of  the  maximum  steady-state  rate  of  flow  that  could  be  sustained 
between  two  specified  vertices,  and  the  determination  of  a  specific  family 
of  minimum-cost  flow  patterns  between  these  vertices — one  pattern  for 
every  integral  rate  of  flow  from  zero  to  the  maximum.  That  problem  is 
essentially  Problem  A  of  this  paper.  The  restriction  on  parallel  links  has 
been  removed  here,  and  the  problem  and  its  solution  are  presented  in 
terms  more  closely  allied  to  the  related  Problems  B  and  C. 

The  method  adopted  in  ORO-TP-15^  was  one  of  discovering  a 
succession  of  path  flows  from  the  origin  to  the  destination  in  such  a  way 
that  each  new  path  flow  contributed  one  more  unit  of  flow  with  the  least 
possible  increase  in  cost.  The  mechanism  for  selecting  the  "best"  path 
at  each  stage  of  the  process  consisted  of  associating  an  economic  length 
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(called  the  ”  effective  length”)  with  each  direction  of  every  link.  This 
length  was  a  function  of  the  appropriate  constant  unit  cost  and  also  of  the 
amount  of  flow  already  present  in  the  link  owing  to  paths  previously 
selected.  By  defining  these  lengths  properly  the  best  path  to  use  next 
was  simply  the  shortest  in  terms  of  total  economic  length. 

A  vital  point  in  the  procedure  was  the  fact  that  a  starting  feasible 
flow  pattern  having  minimum  cost  was  available.  ( This  was  simply  the 
pattern  having  zero  flow  in  every  link  and  consequently  zero  total  cost.  ) 

An  important  extension  of  this  problem  arises  if  one  wishes  to  "  force" 
flow  through  one  or  more  links.  That  is,  the  flow  in  certain  links  is 
required  to  be  in  a  specified  direction  and  to  be  of  a  magnitude  lying 
between  two  positive  bounds.  In  this  case  there  may  not  be  any  feasible 
flow  patterns  if  the  constraints  are  too  stringent.  If  any  feasible  flows 
between  the  specified  origin  and  destination  exist  then  there  are  integers 
m  and  M  that  represent  the  minimum  and  maximum  feasible  rates  of  flow 
between  these  points.  One  can  again  pose  the  problem  of  determining  a 
family  of  minimum-cost  flows,  one  for  each  integral  rate  of  flow  between 
m  and  M.  This  is  Problem  B  of  the  present  paper.  One  can  deal  readily 
with  problems  of  this  type  by  employing  an  appropriate  modification  of 
the  network,  whose  effect  is  to  remove  the  constraints  that  "  force"  flow, 
and  hence  to  reduce  the  problem  to  one  of  the  previous  type.  The  pro¬ 
cedure  adopted  here  involves  a  network  augmentation  that  differs  slightly 
2 

from  that  of  Berge,  together  with  the  necessary  extension  of  ideas 
required  to  find  a  complete  set  of  minimum-cost  flows  rather  than  simply 
one. 
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The  basic  combining  properties  of  elementary  flows  (corresponding 
to  paths  and  cycles),  which  play  an  important  role  in  the  solution  of 
Problem  A  (and,  indirectly,  Problem  B),  also  yield  a  criterion  for 
testing  the  optimality,  in  terms  of  cost,  of  a  given  feasible  flow  pattern. 
Problem  C  of  this  paper  concerns  such  a  test.  A  systematic  method  for 
testing,  and  for  reducing  the  cost  of  nonoptimal  patterns,  is  incorporated 
in  the  resulting  procedure. 

Network  flow  problems  have  been  the  object  of  considerable  attention 

in  the  literature.  The  list  of  references  given  at  the  end  of  this  paper  is 

minimal.  Concepts  and  results  from  these  particular  references  were 

2  3 

used  to  form  the  basis  for  the  present  exposition.  Berge,  Fulkerson, 

4  5 

Ford  and  Fulkerson,  and  Iri  give  procedures  for  finding  minimum-cost 
feasible  flow  patterns.  Of  these  Ford  and  Fulkerson  and  Iri  find  complete 
families  of  minimum-cost  patterns  transmitting  increasing  amounts  of 
flow  from  origin  to  destination.  The  principle  of  developing  "  expanding” 

flow  patterns  by  successive  addition  of  flows  along  appropriate  ” unsaturated" 

6  7 

paths  is  discussed  by  Dantzig  and  Fulkerson.  Pollack  and  Wiebenson 

review  a  variety  of  algorithms  for  finding  shortest  routes  through  net¬ 
works.  The  particular  shortest- route  algorithm  employed  here  (as  the 
repeated  "  inner  loop"  in  the  procedure  for  solving  Problem  A)  is  a 

g 

"labeling"  process  patterned  after  those  presented  by  Ford  and  Ford 
9  10 

and  Fulkerson.  The  cited  publication  by  Jewell  includes  a  development 
of  ideas  similar  to  that  adopted  here,  particularly  for  the  solution  of 
Problem  A.  Since  a  comprehensive  canvass  of  the  literature  of  flow 


5 


theory  was  not  undertaken,  other  related  material  besides  that  referenced 
here  has  appeared  in  the  literature. 

The  next  section  of  this  paper  presents  both  fundamental  terminology 
and  results  required  for  the  precise  formulation  and  solution  of  the  three 
flow  problems  considered.  The  following  section  states  the  three  prob¬ 
lems.  A  basic  vocabulary  having  been  established,  the  next  section 
indicates  how  all  three  mathematical  problems  arise  in  the  contest  of  a 
general  transportation  problem.  (The  reader  who  does  not  wish  to  have 
the  continuity  of  the  mathematical  development  broken  may  omit  this 
section  without  loss  of  any  essential  mathematical  facts.  )  The  remaining 
three  sections  develop  and  justify  the  procedures  for  solving  the  three 
problems.  A  glossary  of  recurring  symbols  and  an  index  of  principal 
terms  are  included  at  the  end. 
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GENERAL  DEFINITIONS  AND  RESULTS 


Let  X  and  U  denote  two  sets,  whose  elements  will  be  referred  to  as 
vertices  and  arcs  respectively.  Assume  that  two  distinct  vertices  xeX 
and  ye  X  are  associated  with  each  arc  ue  U,  and  are  designated  as  the 
initial  and  terminal  vertex  of  u,  respectively.  This  relationship  is 
expressed  symbolically  by  writing  u  =  (x, y).  This  is  read  "arc  u 
extends  from  vertex  x  to  vertex  y,  "  because  of  the  geometric  interpretation 
described  below. 

The  pair  (X,  U)  together  with  the  association  of  vertices  with  arcs  is 
a  directed  graph.  If  X  and  U  are  both  finite  sets  the  graph  is  called  a 
finite  directed  graph.  A  finite  directed  graph  can  always  be  represented 
geometrically  in  three-dimensional  Euclidean  space  by  a  set  of  points 
corresponding  to  the  vertices  and  a  set  of  simple  open  curves*  repre¬ 
senting  the  arcs.  Each  curve  is  constructed  so  that  its  end-points  are 
the  points  in  space  corresponding  to  the  two  vertices  associated  with  the 
arc  it  represents.  Moreover  the  curves  can  be  constructed  in  such  a 

*A  simple  open  curve  is  a  point  set  topologically  equivalent  to  (or 
homeomorphic  witl|  a  straight-line  segment;  i.e. ,  it  can  be  transformed 
into  a  segment  by  a  suitable  one-to-one  transformation  continuous  in  both 
directions. 
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manner  that  they  are  mutually  disjoint,  except  for  common  end-points 
they  may  share.  Finally,  each  curve  is  considered  to  be  directed,  the 
direction  being  induced  by  traversing  the  curve  from  its  initial  to  its 
terminal  vertex.  (Actually,  a  much  larger  class  of  directed  graphs  than 
merely  the  class  of  finite  graphs  has  geometric  realizations  of  the  above 
type.  Any  graph  for  which  neither  X  nor  U  has  a  higher  order  of  infinity 
than  the  set  of  real  numbers  can  be  represented.  )  Most  of  the  terais  and 
concepts  associated  with  graphs  are  particularly  transparent  when  related 
to  geometric  graphs,  as  structures  of  the  above  type  are  called. 

The  initial  and  terminal  vertices  of  an  arc  are  collectively  called  its 
end-points.  An  arc  is  said  to  be  incident  with  its  end-points,  and  non¬ 
incident  with  all  other  vertices.  A  vertex  that  is  not  incident  with  any  arc 
is  an  isolated  vertex.  If  u^  and  u^  are  distinct  arcs  such  that  u^^  s  (x,  y ) 
and  either  U2  =  (x,y)  or  U2  ^  (y,x)  then  Uj^  and  Ug  are  termed  parallel 
arcs.  If  u.,  s  (x,  y)  then  u,  and  u„  are  strictly  parallel  arcs.  Two  dis- 

'  C  1  to 

tinct  arcs  are  adjacent  arcs  if  they  have  a  common  end-point.  They  are 
strictly  adjacent  arcs  if  the  terminal  vertex  of  one  coincides  with  the 
initial  vertex  of  the  other.  If  Uj^  and  are  strictly  adjacent,  and  the 
terminal  vertex  of  u^^  coincides  with  the  initial  vertex  of  u^,  u^^  is  said  to 
precede  Ug. 

A  path  is  a  set  of  arcs  that  can  be  arranged  in  a  sequence 
{  Uj^,  Ug,  .  .  . ,  u^}  such  that  u^  precedes  u.^^^  for  i  =  1,  2,  .  .  . ,  n  -  1 
and  such  that  the  initial  vertex  of  u^^  is  distinct  from  the  terminal  vertex 
of  Uj^.  More  specifically,  the  above  path  is  called  an  x  — y  path  if  x  is  the 
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initial  vertex  of  and  y  is  the  terminal  vertex  of  u^.  In  a  geometric 

directed  graph  an  x  -*y  path  corresponds  to  an  open  curve  with  end-points 

X  and  y  such  that  if  one  envisions  "  traversing"  it  from  x  to  y,  traversing 

the  arcs  in  the  order,  u- ,  u„,  .  .  . ,  u  ,  then  the  direction  of  each  arc 

X  ^  n 

agrees  with  the  orientation  induced  by  this  traversal.  A  path  is  a  simple 
path  if  no  vertex  is  incident  with  more  than  two  arcs  of  the  path.  A  cycle 
is  defined  in  the  same  manner  as  a  path  except  that  the  initial  vertex  of  u^ 
coincides  with  the  terminal  vertex  of  u^.  A  simple  cycle  is  one  such  that 
no  vertex  is  incident  with  more  than  two  arcs  of  the  cycle.  Figure  1 
illustrates  these  important  substructures  of  a  graph. 

It  will  be  convenient  to  consider  also  sets  of  arcs  that  differ  from 
paths  and  cycles  only  in  that  certain  of  the  arcs  have  the  "wrong"  orien¬ 
tation.  If  (X,  U)  is  a  directed  graph  let  U'  denote  a  set  of  elements  that 
are  in  one-to-one  correspondence  with  the  arcs  of  U,  and  let  u'  denote 
the  element  of  U'  corresponding  to  ueU.  The  elements  of  U'  will  also  be 
considered  as  arcs  whose  end-points  are  specified  by  the  relation 
u'  =  (y,x)  whenever  u  =  (x,y).  The  symbol  W  will  always  denote  the 
enlarged  set  U  U  U'.* 

Consider  now  the  directed  graph  (X,  W)  obtained  by  enlarging  (X,  U) 
in  the  above  manner.  If  a  simple  x  -*y  path  {  u^  u^,  •.  .  . ,  Uj^}  is  given 
in  (X,  W),  in  general  some  of  the  arcs  will  be  members  of  U,  and  others 
of  U'.  Arcs  of  the  former  type  are  called  normal  arcs,  relative  to  the 

*A  U  B  denotes  the  union  of  sets  A  and  B;  i.e. ,  the  set  of  elements 
appearing  in  one  or  both  sets. 


9 


o.  Bos  ic  network 


b.  An  X3-»xjt  port) 


d.  A  cycle 


^'9'  l~lllu«trotion  of  Paths  ond  Cycles 


path,  and  arcs  of  the  latter  type  are  called  inverted  arcs.  Similarly, 
given  a  simple  cycle  in  (X,  W)  the  arcs  can  be  classified  as  normal  or 
inverted  in  the  same  way. 

A  directed  graph  (X,  U)  is  said  to  be  connected  if  for  every  pair  of 
distinct  vertices  x  and  y  there  is  at  least  one  simple  x  -•y  path  in  the 
extended  graph  (X,W). 

In  the  remainder  of  this  paper  consideration  will  be  restricted  to 
directed  graphs  of  a  special  type.  Specifically,  a  network  is  defined  as 
a  finite  directed  graph  that  is  connected. 

Let  (X,  U)  be  an  arbitrary  network.  A  flow  pattern,  or  simply  a 
flow  in  (X,  U),  is  a  function  f  that  maps  the  set  of  arcs  U  into  the  set  of 
integers.  If  ueU,  the  integer  f(u)  is  called  the  flow  in  arc  u.  If  u  =  (x, y ) 
the  flow  is  said  to  be  "from  x  to  y,  "  or  "  toward  y  "  and  "  away  from  x,  " 
if  f(u)  >0.  It  is  said  to  be  " from  y  to  x,  "  or  "  toward  x "  and  "away  from 
y,  "  if  f(u)  <0. 

This  terminology  is  suggested  by  interpreting  (X,  U)  as  a  set  of 
essentially  one- dimensional  channels  (arcs)  interconnecting  certain 
locations  (vertices)  where  some  homogeneous  substance  is  flowing 
through  each  arc  u  at  a  steady  rate  f(u)  ,  and  where  the  sign  of  f(u), 
in  conjunction  with  the  orientation  of  u,  indicates  the  direction  of  flow. 

Let  U(-*x)  denote  the  set  of  arcs  of  U  having  x  as  terminal  vertex, 
and  let  U(x-")  denote  the  set  of  arcs  of  U  having  x  as  initial  vertex. 
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Given  a  flow  pattern  f  the  summation*  denotes  the  total 

magnitude  of  flow  away  from  x,  minus  the  total  magnitude  of  flow  toward 
X,  contributed  by  arcs  in  the  set  U(x-^).  Similarly  21y^_^^jf(u)  denotes 
the  total  flow  toward  x,  minus  total  flow  away  from  x,  for  arcs  of  U(-*x). 

The  net  output  at  vertex  x  relative  to  flow  pattern  f.  denoted  by  n^(x),  is 
defined  as  follows: 

If  either  of  the  sets  U(x-*)  or  U(-*-x)  is  void  the  corresponding 
summation  is  taken  to  have  value  zero.  ( It  is  a  consequence  of  the 
connectivity  of  (X,  U)  that  not  both  sets  can  be  void. ) 

One  can  classify  the  vertices  of  (X,  U)  in  terms  of  the  values  of 
their  respective  net  outputs  as  follows: 

If  fij(x)  >0,  x  is  a  source. 

If  flj(x)  <0,  X  is  a  sink. 

If  fl£(x)  =  0,  X  is  an  intermediate  vertex. 

It  is  readily  seen  that  =  0  since,  when  Z;j^ft£(x)  is  expanded, 

the  resulting  double  summation  involves  f(u)  precisely  twice  for  every 
arc  ueU — once  with  a  positive  sign  and  once  with  a  negative  sign.  [  Since 
u€U(x  -•)  for  precisely  one  xeX,  and  ucU(-*y)  for  precisely  one  yeX,  and 
y  X.  ]  Figure  2a  illustrates  the  above  notions. 

If  Cj(x)  =  0,  except  possibly  for  x  =  x^  and  x  =  x^,  then  ^^(x^ )  =  -  fij(Xj ). 
So  either  every  vertex  is  an  intermediate  vertex  or  else  x^  and  x^  are  the 


*Summations  of  the  form  will  be  written  as  simply  2^  when 
ambiguity  will  not  arise. 
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n{*7)-  5 
n(x,)-  2 

0  {*4)  “  -4 
fl  (*j) »  -3 
n  (otKcrt)  -  0 


a.  A  (low  pottorn 


n(xj)-  2 
n(x7)<-2 

n  (othori)  -  0 


n  (x.)  -  1 
n  (xj)  -  -] 

n  (olhort)  -  0 


Fig.  2 — ExampUs  of  Flow  Potlorns 


13 


source  and  sink,  in  some  order.  Assuming  that  is  the  source,  i.e. , 
that  ^^(x. )  >  0,  such  a  pattern  is  called  an  x.  — Xj  flow  pattern.  Given  an 
X.  -•■x^  flow  pattern  f,  fij(x^ )  is  called  the  value  of  the  flow  pattern.  Intu¬ 
itively,  it  represents  the  rate  at  which  flow  "moves"  from  x.  to  x^.  It 
will  be  convenient  to  include,  in  the  class  of  x^  -♦Xj  flow  patterns,  patterns 
f  also  such  that  fi^(x)  =  0  for  all  xcX,  including  x.  and  x^.  Such  patterns 
will  be  considered  as  zero-valued  x^  -♦x^  flow  patterns. 

If  f  and  g  are  two  flow  patterns  defined  on  the  same  network  (X,  U), 
f  and  g  are  said  to  be  conformal  if  there  is  no  arc  ueU  such  that 
f(u)  •  g(u)  <0.  Expressed  differently,  for  every  arc  ueU,  either  f(u)  >  0 

and  g(u)  >  0  or  else  f(u)  <  0  and  g(u)  <  0.  A  set{  f^^,  fg . fj^}  of 

flows  is  conformal  if  the  flows  are  pairwise  conformal. 

The  sum  of  two  flow  patterns  f  and  g,  denoted  by  f  +  g,  is  defined  by 

the  relation  (f  +  g)(u)  =  f(u)  +  g(u)  for  all  uelJ,  Similarly,  the  difference 

f  -  g  is  defined  by  the  relation  (f  -  g)(u)  =  f(u)  -  g(u)  for  all  ueU.  If  f^^ 

and  fg  are  both  x^  -•Xj  flow  patterns  having  values  k^  and  respectively, 

it  is  easily  seen  that  L  +  f„  is  also  an  x.  -•x.  flow  pattern  whose  value  is 

1  ^  1  j 

kj^  +  kg.  Moreover,  fj^  -  fg  is  a  flow  pattern  of  value  |  “  *^2!  ’ 

X.  -"x^  or  x^  -•X.  flow  pattern  if  k^^  -  kg  >  0  or  <  0  respectively. 

Two  particularly  simple  types  of  flow  patterns,  called  elementary 
flow  patterns,  play  a  central  role  in  the  following  development.  It  will 
be  seen  that  these  constitute  basic  building  blocks  for  more  complex  flow 
patterns  in  the  sense  that  a  given  pattern  can  be  decomposed  into  (expressed 
as  the  sum  of)  appropriate  elementary  patterns,  and  patterns  of  the  types 
sou^t  later  can  be  synthesized  from  appropriate  elementary  patterns. 
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Let  (X,  U)  be  an  arbitrary  network,  and  let  (X,  W)  denote  the  enlarged 


network  described  earlier  where  W  =  UUU'.  If  C  is  a  simple  x-»y  path 
in  (X,  W),  an  associated  flow  pattern  f  in  (X,  U)  is  defined  as  follows: 

1  if  ucC 

f(u)  =<-l  if  u’cC 

0  if  neither  ueC  nor  u'cC. 

It  is  easily  verified  that  f  is  an  x-*y  flow  pattern  whose  value  is  1,  Such 
flow  patterns  are  called  elementary  x -^y  path  flows  (see  Fig.  2). 

If  C  is  a  simple  cycle  in  (X,  W),  and  f  is  defined  in  the  same  manner, 
it  is  easily  seen  that  f  is  a  zero-valued  flow  pattern,  called  an  elementary 
cycle  flow. 

The  following  lemma  establishes  the  fact  that  an  arbitrary  x^  -•X2 
flow  pattern  can  be  expressed  as  the  sum  of  conformal  elementary  flows. 


Lemma  1 

If  f  is  an  Xj^  -•Xg  flow  pattern  of  value  k  in  a  network  (X,  U),  and 
f(u)f  0  for  at  least  one  arcucU  then  it  is  possible  to  express  f  as  a  sum 
of  the  form 


f  -  g^  +  .  .  .  +  g^  +  h^  +  .  .  .  +  h^ 

where  the  g^'s  are  k  elementary  Xj^  -^Xg  path  flows,  the  h/s  are 
elementary  cycle  flows,  and  the  g/s  and  h/s  are  jointly  conformal. 

Proof.  Assume  first  that  k  >  0.  Since  )  =  k  >  0  there  exists 

either  an  arc  u€U(x^  -•)  such  that  f(u)  >  0  or  else  an  arc  veU(-*Xj^ )  such 
that  f(v)  <  0.  (Of  course  arcs  of  both  types  may  exist.  )  Let  u^  denote  an 
arc  of  either  type,  and  let  y^^  denote  its  other  end-point.  Since  u^^ 
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contributes  flow  toward  and  since  )  =  0  (assuming  y^^  ^  Xg) 
there  is  either  an  arc  ueU(yj^  -•)  with  f(u)  >  0  or  else  an  arc  v€U(-*yj^> 
with  f(v)  <0.  Let  Ug  be  such  an  arc  and  let  denote  its  other  end¬ 
point.  Continuing  in  this  manner,  one  of  two  possibilities  must  eventually 
arise.  A  stage  k  may  occur  such  that  {  y^^,  y^,  .  .  . ,  yj^  are  all 

distinct  vertices,  but  that  yj^  =  or  =  y^  for  some  j  <  k.  If  this  does 

not  happen  a  stage  k  must  occur  eventually  such  that  y^^  =  Xg,  since  there 

are  a  finite  number  of  vertices.  In  either  case  let{  Wj^,  . Wj^} 

denote  the  sequence  of  arcs  of  W  related  to{  u^^,  u^,  .  .  .,  Uj^}  as 
follows:  w^  =  Uj  if  f(up  >  0  (i.e. ,  if  y^  is  the  terminal  vertex  of  u^ ), 
while  w.  =  u!  if  f(u. )  <  0  (i.e. ,  if  y.  is  the  initial  vertex).  If  y.  =  y.  for 
some  j  <  k  then{  ^^4.2*  •  •  •>  ®  simple  cycle.  Similarly, 

if  yk  ~  x^  entire  set  of  w/s  is  a  simple  cycle.  If  the  y^'s  are  distinct 
but  yj^  =  Xg  then  the  entire  set  of  w^'s  forms  a  simple  x^  path. 

Denote  by  1^^  the  elementary  cycle  flow  or  elementary  x^^  -♦Xg  path 
flow  corresponding  to  the  cycle  or  path  determined  in  this  manner,  and 
set  ®  f  In  the  former  case  f^^  is  again  an  x^^  -►Xg  flow  pattern  of 

value  k.  In  the  latter  case  f^^  is  an  x^^  -•Xg  flow  pattern  of  value  k  -  1. 

In  either  case  if  the  value  of  f^^  is  greater  than  zero  a  second  elementary 

flow  Ig  can  be  "  extracted"  by  repeating  the  process.  Let  ^2  ~  ~  ^  ~ 

i 

1  -  1„.  In  general,  having  determined  f.  =  f  -  2  1.  either  f.  has  value 

1  ^  1  j=l  J  1 

zero  or  else  another  elementary  path  or  cycle  flow  can  be  found  so 

i+1 

that  one  can  define  f.^.,  =  f  -  2  1..  But  every  time  an  elementary  flow 
1+1  1  J 

is  subtracted  the  effect  is  to  reduce  the  value  of  certain  positive  arc 
flows  by  one  unit,  increase  the  value  of  certain  negative  arc  flows  by 
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one  unit,  and  leave  unchanged  the  flows  in  arcs  not  associated  with  the 

cycle  or  path.  Such  modifications  can  occur  only  a  finite  number  of 

times  before  all  arc  flows  are  reduced  to  zero.  So  eventually  a  stage  m 

m 

is  reached  such  that  f  =  f  -2  1.  is  a  zero-valued  flow  pattern  z,  and 

m  1  j 

{  1^,  Ig,  .  .  . ,  1^)  consists  of  k  elementary  -•Xg  path  flows  and  m  -  k 
cycle  flows. 

If  z  is  identically  zero  the  desired  decomposition  has  been  produced. 

If  not,  let  y  be  any  vertex  such  that  z(u)  >0  for  some  ueU(y  -►)  or  else 

z(v)  <  0  for  some  veU(-*y).  Since  fi  (x)  =  0  for  all  xeX,  one  can  again 

z 

generate  a  sequence  of  arcs  as  before.  In  this  case  and  do  not  play 
an  exceptional  role,  and  e'  entually  the  path  determined  by  "tracing" 
flow  through  a  succession  of  vertices  must  lead  back  to  a  vertex  previ¬ 
ously  encountered.  When  this  occurs,  another  cycle  flow  1^+^  is 
determined,  yielding  another  zero-valued  pattern  z^^  =  z  -  Either 

Zj^(u)  =  0  for  all  ueU,  or  still  another  cycle  flow  l^j^+g  "  extracted.  " 


1  ,  where  Z  is 
P 


Ultimately  one  obtains  Z  =  z  -  ~  ^m+2  ’  '  ' 

m  P 

identically  zero.  But  then  f  =  2L  +  z=Z:L  where  k  of  the  L's  are  path 
flows  and  p  -  k  are  cycle  flows.  Setting  n  =  p  -  k  the  desired  decompo¬ 
sition  has  been  obtained.  (If  f  was  initially  of  value  zero,  only  the  part  of 
the  process  related  to  decomposing  z  is  required.  ) 

The  fact  that  the  elementary  flows  are  conformal  follows  from  the 
nature  of  their  selection.  Arc  flows  were  consistently  chosen  that  agreed 
in  direction  with  those  of  the  original  pattern  f.  If  f{ u )  >  0  then  L(u )  >  0 
for  every  j,  while  if  f(u)  <0  then  lj(u)  <  0  for  every  j.  (More  precisely 
in  the  former  case  L(u)  =  1  for  f(u)  of  the  elementary  flows  and  L(u)  =  0 
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for  the  others.  In  the  latter  case  ~  ~1  for  ~f(u)  elementary  flows 
and  lj(  u)  =  0  for  the  rest. )  This  completes  the  proof. 

In  order  to  clarify  the  foregoing  procedure  Fig.  3  illustrates  two 
decompositions  of  a  flow  pattern  of  value  2  into  conformal  elementary 
flows.  (It  should  be  noted  that  such  decompositions  are  not  in  general 
unique.  ) 

A  network  (X,  U)  is  said  to  be  capacitated  if  two  integer-valued 
functions  b  and  c  are  defined  on  U  and  satisfy  the  following  relation: 
c(u)  >  0  and  c(u)  >  b(u)  for  all  ueU  . 

The  integers  c(u)  and  b(u)  are  called  the  upper  and  lower  bounds  on 
flow  in  arc  u  respectively.  A  flow  pattern  f  is  feasible  in  arc  u  if 
b(u)  <  f(u)  <  c(u).  If  b(u)  <  0  this  means  that  the  set  of  feasible  flows 
in  u  range  from  c(u)  units  of  flow  in  the  direction  associated  with  u  to 
-b(  u )  units  of  flow  in  the  reverse  direction.  If  b(  u )  >0  one  says  that  the 
flow  in  u  is  "forced"  to  be  in  the  direction  of  arc  u,  and  to  range  from  a 
minimum  of  b(u)  to  a  maximum  of  c(u)  units  of  flow.  A  flow  pattern  is 
feasible  in  a  network  if  it  is  feasible  in  every  arc  of  the  network. 

The  following  lemma  concerning  feasible  flow  patterns  will  be 
required  in  later  sections. 

Lemma  2 

Let  fj^  and  be  feasible  x-*y  flow  patterns  having  values  k^^  and 
^2-^1’  ^  pattern  fj^  -  fg  be  written  as  a  sum  of 

conformal  elementary  flows: 

■  ^2  ’  *1  ®k  ^ 
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c.  Altarnota  decomposition  t 

Fig.  3— Decomposition  of  o  Flow  Pattern  into  Elements 

tin  b  and  e,  arrows  indicate  direction  of  flow,  not  arc 
oriontotions  as  in  a. 
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where  k  =  -  k^.  Then  +  S  is  also  a  feasible  x  -•y  flow  if  S  is  the 

sum  of  any  subset  of  (  ,  gj^,  h^ . h^}. 

Proof.  Since  the  g.'s  and  K's  are  jointly  conformal,  in  any  arc  u 

they  are  all  nonnegative  if  f^(u)  >  f2(u)  or  all  nonpositive  if  fj^(u)  <  f2(u). 

In  either  case  adding  a  subset  of  them  to  yields  an  arc  flow  f(u)  that 

lies  between  the  values  f,  (u)  and  f„(u).  Since  f,  (u)  and  f„(u)  each  lie 

i.  Z  1  Z 

between  b(u)  and  c(u)  it  follows  that  f(u)  does  also.  Hence  f  is  a  feasible 
flow  pattern.  This  completes  the  proof. 

A  network  is  said  to  be  weighted  if  two  functions  a^  and  a^  called  unit 
cost  functions  are  defined  on  the  set  of  arcs,  whose  values  are  nonnegative 
real  numbers.  If  u  is  any  arc,  aj^(u)  is  considered  the  unit  cost  of  flow 
in  the  direction  of  the  arc,  and  a„(u)  is  the  unit  cost  of  flow  in  the  oppo¬ 


site  direction.  If  f  is  a  flow  pattern  the  cost  of  flow  in  u,  denoted  by 
T  [f(u)],  is  defined  as  aj^(u)  •  f(u)  if  f(u)  >  0,  and  as  -a2(u)  •  f(u)  if 
f(u)  <  0.  Thus  in  any  case  the  cost  of  flow  in  an  arc  is  nonnegative.  The 
total  cost  of  a  flow  pattern  f  is  denoted  by  T(f).  (Two  cost  functions  are 
specified  in  order  that  asymmetric  unit  costs  may  be  reflected.  )  The 
problems  investigated  in  the  remainder  of  this  paper  deal  with  con¬ 
structing  feasible  flow  patterns  whose  costs  are  minimum,  in  a  given 
capacitated  and  weighted  network.  [  Of  course  if  b(u)  >  0,  a2(u)  is 
irrelevant,  since  no  feasible  arc  flow  is  such  that  f(U)  <  0. )  Feasible 
patterns  whose  costs  are  minimum  will  be  called  ideal  patterns. 

An  important  relation  between  the  costs  of  certain  flow  patterns  is 
given  in  the  following  result. 
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Lemma  3 


If  f,  g,  and  h  are  x  -"y  flow  patterns  in  the  same  network,  and  g  and 
h  are  conformal,  then 

T(f  +  g  +  h)  -  T(f  +  g)  >  T(f  +  h)  -  T(f)  . 

Expressed  differently  the  incremental  cost  resulting  from  adding  h  to 
f  +  g  is  greater  than  or  equal  to  the  incremental  cost  resulting  from 
adding  h  to  f  alone. 

Proof.  Let  u  denote  an  arbitrary  arc  of  U,  and  set 

T^  =  T[f(u)  +  g(u)  +  h(u)]  -  T[f(u)+  g(u)]  -  T[f(u)  +  h(u)]  +  T[f(u)]  . 

One  can  also  write  T  as 
u 

Tu  =|f(u)+g(u)+h(u)|a.(u)  -  |f(u)+g(u)|  a.(u)  -  |f(u)+h(u)|  aj^(u)+ |f(u)|  a^(u), 
where  a.(u)  is  a,  (u)  or  a„(u)  depending  on  whether  [  f(u)  +  g(u)  +  h(u)] 
is  >  0  or  <  0,  and  a^,  a^,  and  a^  are  similarly  determined.  It  will  be 
shown  that  >  0  for  all  ueU. 

Suppose  that  f(u)  >  0.  [  If  f(u)  <  0,  u  could  simply  be  reoriented  and 
f(u)  could  be  replaced  by  -f(u)  without  essentially  changing  the  flow  pat¬ 
tern.  Of  course  then  g(u)  and  h(u)  must  also  be  replaced  by  their  nega¬ 
tives  and  aj^(u)  interchanged  with  a^fu).  ]  Since  g  and  h  are  conformal, 
one  of  the  following  sets  of  inequalities  must  hold: 

f(  u )  >  0,  g( u )  >  0,  and  h(  u )  >  0  , 


or 


f(  u )  >  0,  g( u )  <  0,  and  h(  u )  <  0  . 

If  the  first  set  of  inequalities  holds,  then  the  absolute  value  signs  are 
superfluous  in  the  expression  for  T^,  and  a^  =  a^  =  aj^  =  a^  ~ 
evaluating  T^,  =  0  is  obtained.  If  the  second  set  of  inequalities 
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holds,  but  [  f(u)  +  g(u)  +  h(u)]  >  0,  then  [  f(u)  +  g(u)]  >  0  and 
[  f(u)  +  h(u)]  >  0,  and  =  0  is  again  obtained.  The  remaining  case  is 
that  for  which  f(u)  >  0,  g(u)  <  0,  h(u)  <  0,  and  f(u)  +  g(u)  +  h(u)  <  0. 

In  this  case,  can  be  written  as 

Ty  =  -[  f(u)+g(u)+h(u)]  a2(u)  -  |f(u)+g(u)  |  a.(u)  -  |f(u)+h(u)|  aj^(u)  +  f(u)aj^(u). 

Depending  on  the  relative  magnitudes  of  f(u)  and  g(u),  jf(u)  +  g(u)  j  a^(u) 
is  equal  to  either  [  f(u)  +  g(u)]  aj^(u)  or  -(  f(u)  +  g(u)]  a2(u).  Similarly 
jf(u)  +  h(u) I  aj^(u)  can  be  expressed  either  as  [  f(u)  +  h(u)]a^(u)  or  else 
as  -[f(u)  +  h(u)]a2(u).  By  substituting  these  expressions  in  the  above 
equation  for  in  all  (four)  possible  combinations  and  simplifying  the 
resulting  expression,  an  expression  for  that  is  clearly _>0  is  always 
obtained.  Hence  for  every  u€  U,  >  0.  Hence  >  0.  But  this  is 

equivalent  to  the  assertion  of  the  lemma,  and  hence  the  proof  is  complete. 
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STATEMENT  OF  FLOW  PROBLEMS  A,  B,  AND  C 


In  this  section  three  specific  network  flow  problems  are  stated.  These 
are  designated  as  Problems  A,  B,  and  C.  The  remainder  of  the  paper  is 
devoted  to  the  presentation  of  complete  solutions  to  these  problems. 

As  common  assumptions  for  all  three  problems  it  is  assumed  that  a 
capacitated  and  weighted  network  (X,  U)  is  given.  Two  vertices  and 
are  distinguished  from  the  remaining  vertices,  and  are  to  be  considered 
as  the  unique  source  and  unique  sink,  respectively,  of  all  flow  patterns  of 
interest.  In  other  words  the  class  of  flow  patterns  under  consideration  is 
limited  to  the  class  of  x^  patterns.  Retaining  the  symbols  intro¬ 

duced  earlier,  b(u)  and  c(u)  denote  the  lower  and  upper  bounds  of  flow  in 
arc  u,  and  aj^(u)  and  a2(u)  the  unit  costs  appropriate  to  the  two  directions 
that  flow  may  assume  in  arc  u. 

For  Problem  A  it  is  further  assumed  that  b(u)  <  0  for  every  ue  U. 

For  such  a  network,  called  a  Type  A  network,  it  is  evident  that  a  feasible 
-•x^  flow  pattern  having  value  0  exists,  viz. ,  the  pattern  defined  by 
f(u)  =  0  for  all  u€  U.  This  pattern  is  in  fact  ideal,  since  T(f)  =  0  and 
T(g)  >  0  for  every  feasible  flow  pattern  g.  It  follows  from  Lemma  2  that 
if  a  feasible  Xj^  -•x^  flow  pattern  of  value  k  exists  in  a  Type  A  network. 
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then  a  feasible  flow  pattern  exists  having  value  i  for  i  =1,  2,  .  .  . ,  k. 

Since  there  are  at  most  a  finite  number  of  feasible  flow  patterns  because  all 
arc  bounds  are  finite  and  arc  flows  must  be  integers,  this  implies  that  there 
exist  ideal  x^^  -*x^  flow  patterns  having  values  1,  2,  .  .  . ,  k.  The  largest 
k  for  which  feasible  x,  -•x  flow  patterns  exist  will  be  called  the  x.  -►x 
capacity  of  the  network.  In  view  of  these  considerations  the  following 
problem  is  always  meaningful; 

Problem  A.  Given  a  network  (X, U)  that  is  weighted  and  capacitated, 
and  such  that  b(u)  <  0  for  all  uf  U;  develop  a  finite  procedure  that  will 
(a)  determine  the  x^  -►x^  capacity  M  of  the  network,  and  (b)  determine  a 
specific  ideal  x^^  -►x^  flow  pattern  having  value  i  for  i  =  1,  2,  .  .  . ,  M. 

By  a  "finite  procedure"  is  meant  one  that,  for  every  network  of  the 
type  under  consideration,  will  yield  a  solution  after  a  finite  number  of 
"simple"  steps.  A  simple  step  is  an  elementary  arithmetical  operation 
or  comparison  of  magnitudes  performable  by  a  digital  computer.  Although 
the  procedure  given  here  will  not  be  broken  into  ultimate  operations  it 
will  be  given  in  sufficient  detail  so  that  this  could  readily  be  done. 

Problem  A  is  discussed  and  solved  in  full  detail  in  ORO-TP-15.^  A 
restatement  of  the  solution  procedure  and  some  of  its  important  charac¬ 
teristics  is  given  in  this  paper  because  of  the  intimate  relation  of  this 
problem  with  Problems  B  and  C. 

The  assumptions  for  Problem  B  differ  from  those  of  Problem  A  in 
only  one  respect:  it  is  assumed  that  b(u)  >  0  for  at  least  one  arc  ue  U. 
Networks  of  this  type  are  called  Type  B  networks.  Thus  there  is  at  least 
one  arc  in  which  the  flow  is  required  to  be  in  a  specified  direction,  viz. , 
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the  direction  associated  with  the  arc.  It  is  easy  to  construct  networks  of 
this  type  for  which  no  feasible  flow,  of  any  value,  exists.  Suppose,  for 
example,  that  U(-*x^)  and  U(x^-*)  are  empty  sets,  so  that  no  arc  is  directed 
toward  the  source  or  away  from  the  sink.  Then  if  ^ 

^u€U(-*x  )  feasible  x^  -"x^  flow  pattern  is  possible  because  the 

first  summation  is  a  lower  bound  on  the  value  of  flow,  whereas  the  second 
summation  is  an  upper  bound.  However,  if  any  feasible  x^ 
patterns  exist,  then  it  will  be  shown  that  there  exist  minimum  and  maxi¬ 
mum  values,  m  and  M,  such  that  feasible  flows  exists  for  every  value  i 
such  that  m  <  i  <  M.  The  following  problem  is  therefore  always  meaningful: 

Problem  B.  Given  a  network  (X,  U)  that  is  weighted  and  capacitated, 
and  such  that  b(u)  >  0  for  at  least  one  arc  ue  U;  develop  a  finite  procedure 
that  will  (a)  determine  whether  any  feasible  Xj^  -►x^  flow  patterns  exist, 

(b)  determine  m  and  M  if  any  feasible  patterns  exist,  and  (c)  determine  a 
specific  ideal  x^^  -•x^  flow  pattern  having  value  i  for  i  =  m,  m  + 1,  .  .  . ,  M. 

Problem  C  differs  from  the  preceding  ones  in  that  it  is  concerned  with 
ideal  and  nonideal  x^^  -*x^  flow  patterns  having  the  same  value.  Since  no 
restrictions  are  placed  on  the  sign  of  b(u),  as  in  Problems  A  and  B,  the 
problem  is  posed  for  both  Type  A  and  Type  B  networks. 

Problem  C.  Given  a  network  (X,  U)  that  is  weighted  and  capacitated, 
and  given  a  feasible  x^  -•x^  flow  pattern  f^  of  value  k;  develop  a  finite 
procedure  that  will  (a)  determine  whether  or  not  f^  is  ideal,  and  (b)  if  f^ 
is  not  ideal,  find  an  appropriate  succession  of  transformations,  starting 
with  fg,  each  of  which  reduces  the  cost  of  the  preceding  flow  pattern  while 
retaining  its  value  and  feasibility. 
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A  final  note  on  the  three  problems  is  appropriate.  A  reasonable 
solution  to  each  problem  must  consist  of  something  substantially  more 
efficient  than  an  exhaustive  procedure  for  generating  and  comparing  all 
possible  feasible  flow  patterns.  Such  a  procedure  would  in  fact  be  finite, 
but  it  is  intended  here  that  a  procedure  be  devised  that  may  be  reasonably 
implemented. 
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A  TRANSPORTATION  APPLICATION  EMBODYING 
PROBLEMS  A,  B,  AND  C 


This  section  considers  a  general  class  of  problems  arising  in  the 
field  of  transportation — problems  having  to  do  with  the  determination  of 
"optimal"  steady-state  flows  of  materiel  between  certain  origins  and 
destinations  interconnected  by  a  transportation  network.  Problem  A  is 
presented  as  the  mathematical  formulation  of  the  basic  problem  to  be 
solved.  Problem  B  then  arises  in  connection  with  the  attempt  to  gain 
control  of  the  "continuity"  of  a  succession  of  related  problems  of  Type  A. 
Finally,  problems  of  Type  C  arise  when  one  wishes  to  alter  some  of  the 
parameters  of  the  problem  and  determine  whether  or  not  optimality  is 
disturbed,  and  (if  it  is )  to  restore  optimality  in  a  computationally 
efficient  manner.  The  mechanics  of  solving  Problems  A,  B,  and  C  will 
not  enter  into  the  discussion.  It  is  simply  assumed  here  that  efficient 
algorithms  are  available. 

For  the  sake  of  concreteness  suppose  that  one  is  concerned  with  the 
problem  of  moving  goods  from  three  locations  y^^,  y^,  and  y^  to  two 
locations  z,  and  z„  through  a  highway  network  whose  configuration  is 
shown  in  Fig,  4a^  (Although  a  very  simple  situation  is  depicted  here,  it 
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Fig.  4 — Schematic  of  a  Transportation  Network 
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should  be  noted  that  in  practice  the  relevant  network  may  have  thousands 
of  arcs.  It  is  in  such  situations  that  the  need  for  powerful  general 
algorithms  becomes  apparent, ) 

Suppose  that  destinations  and  z^  must  receive  and  Rg  tons/day 
respectively,  and  that  origins  y^^,  y^,  and  y^  can  send  at  most  S^,  S^, 
and  Sg  tons/day  respectively.  It  is  assumed  that  S 
so  that  it  is  feasible  to  support  the  requirements  at  the  destinations  unless 
bottlenecks  in  the  highway  system  prohibit  sufficient  flow. 

As  a  constraint  on  the  highway  system  itself,  assume  that  each  link 
joining  two  successive  junction  points  has  associated  with  it  an  integer 
representing  the  maximum  rate  of  flow,  in  tons  per  day,  that  ;an  be 
sustained  during  the  time  interval  being  analyzed.  It  is  assumed  that 
these  link  "capacities"  apply  in  either  direction  and  that  they  allow  for 
retrograde  movement  of  vehicles.  ( For  the  application  considered  here, 
only  the  "forward"  movement  of  goods  from  origins  to  destinations  is 
considered  explicitly.  Retrograde  movement  is  assumed  to  retrace  the 
forward  routes.  ) 

As  a  criterion  for  optimality,  assume  that  it  is  desirable  to  determine 
a  movement  pattern  which  minimizes  total  ton-mileage  per  day.  Thus  the 
"unit  cost"  associated  with  each  link  is  simply  its  physical  length 
measured  in  miles,  so  that  the  "cost"  of  flow  in  a  link  has  the 
dimensionality  of  ton-miles  per  day. 

The  problem  is  then  one  of  determining  the  direction  and  rate  of 
flow  in  each  link  so  that  (a)  the  net  input  to  destination  is  precisely 
R.,  (b)  the  net  output  from  origin  y^  is  at  most  S.,  (c)  the  net  output  is 


i  +  Sg+Sg^R^  +  Rg, 
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zero  at  all  intermediate  junction  points,  (d)  the  flow  in  each  link  does 
not  exceed  its  specified  capacity,  and  (e)  for  all  possible  flow  systems 
satisfying  a  to  d  the  one  selected  involves  minimum  total  ton-mileage  of 
movement  per  day. 

This  problem  can  be  reformulated  as  one  involving  a  single  origin 
(source)  and  destination  (sink)  by  proper  augmentation  of  the  highway 
system  to  include  certain  hypothetical  links.  Figure  4b  shows  the  con¬ 
figuration  of  the  augmented  system.  A  hypothetical  source  is  added, 
together  with  arcs  directed  from  to  y^^,  y^,  and  yg  respectively.  These 
are  assigned  zero  unit  cost  (or  length)  and  capacities  equal  to  the  corre¬ 
sponding  capabilities  ,  S„,  and  S_.  Similarly  a  hypothetical  sink  x  is 
added,  with  arcs  directed  from  z^  and  z^  to  x^.  These  also  are  assigned 
zero  length  and  capacities  equal  to  and  Rg  respectively.  Note  that  the 
links  of  the  original  highway  system  have  also  been  oriented,  so  that  this 
finite,  connected  system  of  vertices  and  arcs  forms  an  network 

in  the  technical  sense.  The  orientations  here  are  for  frame-of-reference 
purposes,  so  that  the  direction  of  flow  may  be  indicated  by  the  algebraic 
sig^.  (Although  the  direction  of  flow  can  usually  be  predicted  beforehand 
for  most  arcs,  in  some  of  the  lateral  arcs  it  is  not  always  clear  which 
direction  an  optimal  pattern  will  require.  ) 

For  arcs  of  the  original  network,  set  b(u)  =  -K  and  c(u)  =  K,  where 
K  is  the  capacity  of  the  corresponding  link,  and  set  aj^(u)  =  ^2(0)  =  L, 
where  L  is  the  length  of  the  link.  For  arcs  of  the  form  v  =  set 

b(v)  =  0,  c(v)  =  S.,  a-(v)  =  0,  and  a„(v)  =  0.  For  arcs  of  the  form 

lx  ^ 

w  s  (z.,x  )  set  b(w)  =  0,  c(w)  =  R.,  a  (w)  =  0,  and  a„(w)  =  0. 

]  n  )  i  z 
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Since  b,  c,  a, ,  and  a„  have  been  defined  for  the  entire  x,  -•x  network, 
12  In’ 

and  the  network  is  of  Type  A  ( since  b(u)  <  0  for  all  u),  Problem  A  can  be 
posed  for  this  network.  The  -►x^  capacity  M  is  clearly  such  that 
M  <  +  Rg,  since  flow  can  enter  x^  only  through  the  two  arcs  that  join 

Zj  and  Zg  to  x^,  and  these  have  capacities  R^  and  Rg  respectively.  Hope¬ 
fully  the  solution  of  Problem  A  will  reveal  that  M  =  R^^  +  Rg,  so  that  the 
requirements  can  in  fact  be  feasibly  fulfilled.  Even  if  they  cannot  the 
maximum- valued  ideal  flow  pattern  will  indicate  the  maximum  tonnage 
that  can  be  feasibly  transmitted  per  day,  and  will  further  indicate  a 
specific  pattern  of  arc  flows  that  realize  this  maximum  tonnage  with 
minimum  cost  (ton-mileage  per  day).  The  flows  in  a  hypothetical  arc 
should,  of  course,  be  interpreted  as  the  amount  of  materiel  sent  from  an 
origin  or  received  by  a  destination,  whichever  is  appropriate.  Intro¬ 
duction  of  these  arcs  is  merely  an  artifice  that  enables  one  to  treat  y^^, 

^2’  ^3’  ^1’  ^2  intermediate  vertices  for  which  input  equals  out¬ 

put,  rather  than  as  five  exceptional  vertices  with  distinctive  constraints. 
Since  the  artificial  arcs  were  assigned  zero  unit  costs  the  cost  of  an 
ideal  flow  pattern  in  the  augmented  network  is  the  same  as  the  cost  of 
the  "  real"  flow,  viz. ,  the  flow  pattern  obtained  by  finally  disregarding 
the  artificial  arcs. 

Suppose  now  that  one  is  faced  with  the  problem  of  determining 
optimum  routings  for  several  successive  time  intervals.  For  example, 
suppose  that  during  month  1  the  average  daily  demands  at  Zj^  and  are 
R^  and  R2,  and  that  during  month  2  they  change  to  Rj^  and  R^  respectively. 
One  can  consider  each  month  independently  as  a  problem  of  Type  A  in  the 
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augmented  network  described  above.  However,  it  may  well  happen  that 
the  optimum  flow  patterns  determined  for  the  two  months  display  a  serious 
lack  of  "  continuity.  "  For  instance,  one  pattern  may  rely  heavily  on 
northern  routes  and  the  other  on  southern  routes.  Or  one  may  rely 
mostly  on  y^^  and  as  sources  of  supply  and  the  other  largely  or  perhaps 
exclusively  on  y^  and  y^. 

Assume  that  a  specific  arc  is  relied  on  heavily  in  the  optimal  solution 
obtained  for  month  1,  and  that  the  arc  is  not  utilized  in  month  2.  Although 
the  solutions  individually  satisfy  the  criterion  of  minimum  ton-mileage  per 
day,  taken  together  they  may  yield  an  unsound  operational  solution  in 
the  light  of  other  considerations  not  reflected  in  the  mathematical  model. 
For  instance,  the  use  of  this  arc  in  month  1  may  imply  considerable 
investment  in  terms  of  road  or  bridge  rehabilitation  if  the  problem 
arises  in  a  military  context.  To  rehabilitate  the  link  and  then  abandon 
it  in  a  later  month  may  well  be  unreasonable  in  view  of  the  overall 
problem. 

As  a  first  measure  month  1  might  be  reanalyzed,  leaving  all  factors 
unchanged  except  for  the  elimination  of  the  arc  in  question.  Suppose 
that  this  is  done,  and  that  the  new  solution  is  not  acceptable,  either 
because  the  minimum  ton-mileage  per  day  is  unacceptably  higher  than 
it  was  before  or  else  because  removal  of  the  arc  creates  a  bottleneck 
preventing  satisfaction  of  the  total  demands  at  and  Zg.  In  any  event 
assume  that  an  analysis  of  the  situation  leads  to  the  conclusion  that  this 
highway  link  must  be  rehabilitated  and  used  in  month  1,  and  that  it  is 
highly  desirable  to  continue  its  utilization  in  month  2,  thus  lessening 
reliance  on  other  portions  of  the  highway  system. 
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One  way  to  handle  this  problem  is  to  reanalyze  month  2,  forcing  flow 
through  the  relevant  arc.  In  other  words  a  Type  B  problem  can  be  posed 
that  differs  from  the  original  Type  A  analysis  of  month  2  only  in  that 
b(u)  >  0  for  this  arc  and  a  certain  minimum  flow  is  required.  Assume 
that  this  analysis  is  carried  out,  obtaining  a  flow  pattern  f  with  cost 
(ton-mileage  per  day)  T(f'),  which  can  be  compared  with  the  original 
month  2  solution  f  having  cost  T(f). 

Since  the  month  2  problem  was  modified  only  in  the  direction  of 
adding  additional  constraints,  clearly  T(f')  >  T(f).  It  can  happen  that 
the  new  solution,  forced  to  involve  the  arc  in  question,  is  as  cheap,  in 
terms  of  ton-mileage  per  day,  as  the  old  one.  This  can  occur  because 
the  ideal  flow  pattern  having  a  given  value  is  not  necessarily  unique,  so 
that  although  the  original  analysis  of  month  2  yielded  an  ideal  pattern  that 
did  not  involve  the  arc,  another  ideal  pattern  may  actually  exist  that 
does.  In  any  case  a  good  "hard"  number  T(f')  -  T{f)  is  obtained,  which 
is  a  measure  of  the  additional  cost  incurred  if  preserving  continuity  of 
flow  in  the  arc  is  insisted  on.  This  is  just  one  example  of  how  a  Type  B 
problem  can  be  posed  to  "  gain  control"  of  a  set  of  solutions  to  Type  A 
problems  when  the  initial  solutions  are  not  sufficiently  compatible  in 
terms  of  important  constraints  that  were  not  initially  built  into  the 
formulation. 

It  is  sometimes  also  profitable  to  formulate  problems  of  Type  C  in 
the  context  of  a  general  transportation  application  of  the  above  t3q)e. 
Suppose  that  there  is  serious  doubt  concerning  the  accuracy  of  some  of 
the  data.  As  an  example,  suppose  that  there  is  good  reason  to  question 
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the  assumed  capacity  of  an  arc,  and  that  this  arc  appears  importantly  in 
the  solution  as  it  now  stands.  Letting  c(u)  be  the  capacity  used  in  the 
analysis,  suppose  that  the  optimal  pattern  saturates  this  arc,  so  that 
f(u)  =  c(u).  Assume  further  that  there  is  now  reason  to  believe  that  2c(u) 
is  a  better  estimate  of  the  arc's  capacity,  either  because  the  original 
estimate  was  bad  or  else  because  appropriate  physical  measures  can  in 
fact  double  the  capacity  without  undue  expense.  One  approach  would  be 
to  start  from  scratch,  posing  a  new  Problem  A  that  differs  from  the 
original  analysis  only  in  that  the  capacity  of  u  is  increased.  Another 
possibility,  and  one  that  may  save  a  great  deal  of  effort,  arises  by 
taking  the  original  solution  f,  known  to  be  optimal  before,  and  inquiring 
whether  or  not  it  remains  optimal  when  c(u)  is  doubled.  Clearly  f  is 
still  a  feasible  pattern,  which  may  or  may  not  be  optimal.  Studying  this 
Type  C  question,  either  f  will  emerge  optimal  in  the  modified  problem, 
or  else  f  can  be  used  as  the  starting  point  for  finding  an  ideal  pattern  f. 
Here  again  a  quantitative  measure  T(f)  -  T(f')  is  obtained  of  the  difference 
it  makes  whether  c(u)  or  2c(u)  is  the  capacity  of  u. 

The  preceding  remarks  are  only  intended  to  be  suggestive  of  a  host 
of  problems  that  can  be  studied  by  using  formulations  of  Types  A,  B,  and 
C  singly  or  in  combinations.  Network  flow  theory  is  applicable  to  many 
problems,  including  problems  in  which  the  network  is  itself  entirely  an 
artifice  and  has  no  physical  counterpart.  Where  the  network  does  exist, 
as  in  transportation  problems,  it  is  especially  easy  to  formulate  the 
appropriate  mathematical  questions  because  the  model  is  so  direct  a 
reflection  of  reality.  If  the  resulting  questions  are  of  any  of  these  three 
types  the  methods  of  this  paper  or  comparable  methods  are  applicable. 
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SOLUTION  OF  PROBLEM  A 


In  this  section  a  resume  of  the  procedure  for  solving  Problem  A  that 
appeared  earlier  in  ORO-TP-15^  is  given,  together  with  the  intermediate 
mathematical  results  required  to  justify  the  procedure.  This  information 
is  repeated  here  because  of  the  intimate  relation  between  Problems  A,  B, 
and  C.  It  will  be  seen  that  the  solution  to  Problem  B  involves  solution  of 
Problem  A  in  a  related  network  and  that  the  solution  of  Problem  C  requires 
application  of  the  central  ideas  that  were  needed  for  solving  Problem  A. 

Recall  that  a  Type  A  network  is  a  weighted  and  capacitated  network 
(X,  U)  such  that  b(u)  <  0  for  every  arc  ue  U, 

It  was  noted  earlier  that  a  Type  A  network  has  the  property  that  there 
exists  an  integer  M,  called  its  -•x^  capacity,  such  that  feasible  x^^  -►x^ 
flow  patterns  exist  having  values  0,  1,  2,  .  .  . ,  M,  but  none  exist  having 
values  greater  than  M. 

Let  (X,  W)  be  the  enlarged  network  obtained  by  introducing  the  set  of 
"reverse"  arcs  U'  in  the  manner  described  earlier.  If  f  is  a  feasible 
-*x^  flow  pattern  in  (X,  U)  and  P  is  a  simple  x^^  -►x^  path  or  a  simple 
cycle  in  (X,W),  P  is  said  to  be  saturated  relative  to  f  if  it  contains  an 
arc  u£  U  such  that  f(u)  =  c(u)  or  an  arc  u’e  U'  such  that  f(u)  =  b(u). 
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otherwise,  P  is  considered  unsaturated  relative  to  f.  If  gp  is  the 
elementary  path  flow  or  cycle  flow  corresponding  to  P  then  it 

follows  that  f  +  gp  is  feasible  if  and  only  if  P  is  unsaturated  relative  to  f. 

The  following  result  provides  the  basis  for  a  general  procedure  for 
finding  feasible  patterns  having  values  1,  2,  .  .  M. 


Lemma  4 

If  f,  is  a  feasible  x^^  -►x^  flow  pattern  in  (X,  U)  having  value  i,  and 
i  <  M  where  M  is  the  x^^  network  then  there  exists  an 

unsaturated  simple  x^^  -^x^  path  in  (X,W).  If  i  =  M  every  simple  -^x^ 
path  is  saturated. 


Proof.  The  last  statement  follows  immediately.  For  assume  i  =  M 
and  P  is  an  unsaturated  simple  x^^  -•x^  path.  Then  f  +  gp  is  feasible, 
where  gp  is  the  corresponding  elementary  x^^  -►x^  flow.  But  ^  +  gp  ^‘as 
value  i  +  1  =  M  +  1,  contradicting  the  assumption  that  M  is  the  x^^  -►x^ 
capacity  of  ( X,  U ) , 


Assume  now  that  0  <  i  <  M,  and  that  f.  is  a  feasible  x^  -♦x^  flow  pattern 
of  value  i.  Let  be  any  feasible  x^  -*x^  flow  pattern  of  value  i+1.  Then 
f -  f .  is  an  -•x^  flow  pattern  of  value  1 .  According  to  Lemma  1 


-  L  =  g^  +  h^  +  hg  + 


.  +  h 


where  g^^  is  an  elementary  x^^  -►x^  path  flow.  But  from  Lemma  2  it  follows 
that  L  +  gj^  is  feasible.  Hence  g^^  must  correspond  to  an  unsaturated  path 
in  (X,  W).  This  completes  the  proof. 

Thus  if  one  has  an  algorithm  for  finding  an  unsaturated  simple  x^  -♦x^ 
path  relative  to  a  given  feasible  -•x^  flow  pattern,  whenever  one  exists, 


36 


the  following  procedure  will  generate  feasible  patterns  for  all  possible 
values: 

(a)  Let  fp  be  the  feasible  flow  pattern  of  value  zero  defined  by 

fQ(u)  =  0  for  all  u€  U. 

(b)  In  general,  having  determined  a  feasible  pattern  L  of  value  i,  apply 
the  algorithm  to  find  (if  possible)  an  unsaturated  x^  -►x^  path  P. 

(c)  If  no  such  path  P  exists,  i  is  the  x^  capacity  of  the  network, 

and  {  f^,  f^,  .  .  . ,  L  )  is  a  complete  set  of  feasible,  -•x^  patterns. 

(C)  IfP  is  an  unsaturated  x^  path,  let  gp  be  the  corresponding 

elementary  path  flow.  Then  =  L  +  gp  is  a  feasible  pattern  of  value 
i  +  1,  and  step  (b)  can  be  repeated  with  in  place  of  L. 

The  solution  of  Problem  A  revolves  around  the  fact  that  if  care  is 
taken  in  the  above  procedure  to  always  select  the  "best"  unsaturated  path 
at  each  stage  then  the  sequence  (  f^,  f^^,  .  .  , ,  f^^)  will  consist  of  ideal 
flow  patterns.  The  following  property  possessed  by  ideal  flow  patterns 
can  be  established. 


Lemma  5 

If  f.  is  an  ideal  x,  -•x  flow  pattern  of  value  i,  and  i  <  M,  then  an 
elementary  x^^  -•x^  path  flow  g  exists  such  that  f.^j^  =  f.  +  g  is  an  ideal 
Xj^  -*x^  flow  pattern  of  value  i  +  1. 

Proof.  Let  f.^^  be  any  ideal  x^  -►x^  flow  pattern  of  value  i  +  1.  Such 
a  pattern  must  exist  since  i  <  M.  Then  f.^j^  -  L  is  an  x^  -•x^  flow  pattern, 
which  can  ( using  Lemma  1 )  be  written  as 


f. 


i+1 


-  f .  =  +  h,  +  .  .  .  +  h  , 

1  1  n  ’ 
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where  gp  is  an  elementary  path  flow  corresponding  to  some  simple 

path  P,  and  gp  is  conformal  with  the  elementary  cycle  flows  h^^,  hg,  .  .  . ,  h^. 

For  convenience  set  h  =  Sh.,  so  that 

i’ 

f.. ,  =  f.  +  g„  +  h  , 

1+1  1  • 

where  h  is  now  a  zero-valued  flow  pattern  conformal  with  gp.  From 
Lemma  2  it  is  known  that  L  +  gp  and  f ^  k  are  feasible  flow  patterns 
whose  values  are  i  +  1  and  i  respectively.  Lemma  3  can  be  applied  to 
obtain 

T(f.  +.  gp  +  h)  -  T(f.  +  gp)  >  T(f.  +  h)  -  T(fp  . 

Since  f.  is  ideal  the  right-hand  side  is  nonnegative.  Hence 

T(fi  +  gp  +  h)  >  T(f.  +  gp)  . 

But  f.  +  gp  +  h  =  Since  is  also  ideal,  equality  must  hold  in  the 

preceding  expression.  So  f.  +  gp  also  minimizes  T,  i.e. ,  ^ j  gp  ideal 
for  value  i  +  1.  This  completes  the  proof. 

Lemma  5  is  of  fundamental  importance  for  the  attack  on  Problems  A, 

B,  and  C  adopted  here,  and  hence  some  amplifying  remarks  concerning  its 
significance  are  in  order.  If  f.  is  ideal  but  not  maximal,  i.e. ,  if  i  <  M, 
all  possible  feasible  patterns  of  value  i  +  1  can  be  obtained  from  f.  by 
adding  appropriate  patterns  of  value  1.  Specifically,  if  f.^j^  has  value 
i  +  1,  one  merely  adds  -  L)  to  f..  Although  there  are  a  finite  number 
of  distinct  feasible  patterns  of  value  i  +  1,  for  nontrivial  problems  in 
complex  networks  the  number  may  be  astronomical.  Lemma  5  asserts 
that  of  all  patterns  of  value  1  it  is  sufficient  to  consider  only  the  very 
special  class  of  elementary  Xj^  -•x^  path  flows.  Of  these,  only  those 
corresponding  to  unsaturated  paths  are  relevant  since, if  P  is  saturated. 
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relative  to  f.,  then  f.  +  g-  is  not  feasible  and  hence  not  ideal.  Of  course 
the  set  of  unsaturated  paths  may  still  be  extremely  large.  Fortunately  a 
fairly  efficient  constructive  process  can  be  described  for  determining  one 
that  minimizes  T(F  +  gp). 

The  determination  of  the  "best"  -•x^  path  consists  essentially  of 
evaluating  for  each  arc  w  of  U  and  U'  the  Increment  of  cost  that  would 
result  if  one  unit  of  flow  were  added  to  f(w)  and  of  selecting  a  path  of  arcs 
whose  total  incremental  cost  is  minimum.  Before  making  this  precise  it 
is  necessary  to  introduce  several  new  terms. 

If  u€  U,  and  u  =  (x,  y)  let  u  denote  the  elementary  x  -►y  flow  pattern 
obtained  by  considering  u  as  an  x  -^y  path.  Thus  for  all  vc  U 


A 

u(v) 


fl  if  V  =  u 
'^0  if  V  ^  u 


Similarly  let  u*  denote  the  elementary  y  -•x  flow  pattern  corresponding  to 
the  y  -♦x  path  u*,  so  that  for  all  vc  U 


^1/  X  if  v'  ®  u' 
“  «  0  If  V'  f  u' 


Thus  with  every  uc  U  and  every  u’  e  U*  is  associated  the  flow  pattern  u  or 
u',  called  the  characteristic  flow  of  arc  u  or  u'.  If  P  is  a  simple  x^^  -*x^ 
path  the  corresponding  elementary  flow  pattern  gp  can  be  written  as 


W€P 


A 

W. 


If  f  is  a  feasible  x^^  fi®"'  pattern,  and  w  c  U  or  w  €  U’,  the  effective 

cost  of  w.  relative  to  f.  is  denoted  by  e(w;f)  and  defined  as  follows: 


e(w 


-  _  f  T(f  +  w)  -  T(f)  if  (f  +  w)  is  feasible 
’  '  "1  «>  if  (f  +  w)  is  infeasible 
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If  ws  (x, y),  e(w;f)  may  be  considered  as  the  increase  in  cost  which  results 
by  taking  a  feasible  pattern  f  and  superimposing  one  unit  of  flow  in  w  from 
X  to  y.  In  case  this  results  in  an  infeasible  pattern  f  +  w,  the  convention 
is  adopted  that  an  infinite  increase  in  cost  results. 

If  one  considers  the  definition  of  the  cost  of  flow  in  an  arc  in  terms  of 
aj^  and  ag  then  e(w;f)  may  be  evaluated  as  follows: 

For  U6  U, 


e(u;f) 


if  f(u)  =  c(u) 

(u)  if  0  <  f{u)  <  c(u) 
(u)  if  f(u)  <  0 


For  u'€  U', 


e(u';f) 


iff(u)  =  b(u) 

.  a2(u)  if  b(u)  <  f(u)  <  0 
-a^(u)  if  0  <  f(u) 


If  P  is  a  simple  x^  -►x^  path  or  a  simple  cycle  in  (X,  W)  the  effective 
cost  of  P.  relative  to  a  feasible  x^^  — x^  flow  pattern  f,  is  defined  by 


e(P;f)  =  e(w;f)  . 
One  can  establish  the  following  result: 


Lemma  6 

If  f  is  a  feasible  x^^  — x^  flow  pattern  in  (X,  U),  P  is  a  simple  x^  — x^ 
path  or  a  simple  cycle  in  (X,W),  and  gp  is  the  corresponding  elementary 
path  flow, then  e(P;f)  =  «>  if  and  only  if  P  is  saturated,  i.e. ,  if  and  only  if 
f  +  gp  is  infeasible.  If  e(P;f)  <  «  so  that  f  +  gp  is  feasible  then 

T(f  +  gp)  =  T(f)  +  e(P;f)  . 
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Proof.  If  e(P;f)  =  «o  then  e(w;f)  =  «>  for  some  wcP.  Ifw  =  ueU  then 
f(u)  =  c(u).  If  w  =  u'  e  U'  then  f(u)  =  b(u).  In  the  first  case,  u  is  saturated 
in  the  direction  of  its  orientation.  In  the  second  case  it  is  saturated  in  the 
reverse  direction.  In  either  case  it  is  saturated  in  the  direction  associated 
with  P,  so  P  is  saturated.  Conversely  if  P  is  saturated  there  is  an  arc 
w€  P  whose  flow  cannot  be  feasibly  increased  by  one  unit  in  the  direction 
of  w.  If  w  =  u€  U  this  implies  that  f(u)  =  c(u).  If  w  =  u'  e  U'  it  implies 
that  f(u)  =  b(u).  In  either  case  e(w;f)  =  <»;  hence  e(P;f)  =  <». 

It  remains  to  show  that  if  e(P;f)  <  «>  then  T(f  +  gp)  =  T(f)  +  e(P;f). 

Let  w  be  an  arc  of  P,  and  assume  first  that  w  =  ue  U.  Let  T^(f)  denote 
the  cost  of  flow  in  u  relative  to  flow  pattern  f,  and  similarly  for  T^(f+gp). 

If  f(u)  >  0  then  (f  +  gp)(u)  =  f(u)  +  1,  Ty(f)  =  f(u)a^,  and  T^(f+  gp)  =  [  f(u)  + 
l]a^(u).  Hence  T^(f  +  gp)  -  T^(f)  =  aj^(u)  =  e(u;f).  On  the  other  hand  if 
f(u)  <  0  then  (f  +  gp)(u)  =  f(u)  +  1,  =  -f(u)a2(u),  and  T^Cf  +  gp)  = 

-(  f(u)  +  l]a2(u).  Hence  T^^(f  +  gp)  -  Ty(f)  =  -  a2(u)  =  e(u;f).  The 
reasoning  is  similar  if  w  =  u'cU'. 

Hence  for  every  u€  P,  Ty(f  +  gp)  =  T^(f)  +  e(u;f),  and  for  every  u'  €  P, 
T^(f  +  gp)  =  T^(f)  +  e(u';f).  Adding  the  right-hand  side  of  all  these  equations 
yields  T'(f)  +  e(P;f),  where  T'(f)  is  the  total  cost  of  flow  relative  to  f  in 
all  arcs  involved  in  P.  Adding  the  left-hand  members  yields  T'(f  +  gp), 
where  this  e3q)ression  has  a  similar  meaning  but  relative  to  f  +  gp.  Since 
the  cost  of  flow  in  arcs  not  associated  with  P  is  the  same  relative  to  f  and 
to  f  +  gp  it  follows  that  T(f  +  gp)  =  T(f)  +  e(P;f).  This  completes  the 
proof. 
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It  was  seen  earlier  (Lemma  5)  that  if  f  is  a  nonmaximal  ideal  flow 
pattern  then  there  exists  a  path  P  such  that  f  +  gp  is  ideal.  In  view  of 
Lemma  6  one  must  clearly  select  P  in  such  a  way  that  e(P;f)  is  minimum. 

Assume  for  the  moment  that  an  algorithm  can  be  devised,  which  will 
be  called  Algorithm  1,  that  will  perform  the  following  function;  given  an 
ideal  -•x^  flow  pattern  L  of  value  i,  it  will  find  a  simple  x^  -♦x^  path 
P  such  that  e(P;f)  is  minimum  with  respect  to  all  such  paths. 

In  terms  of  Algorithm  1,  which  will  be  described  in  detail  shortly, 
one  can  state  the  following  procedure  that  will  always  solve  Problem  A. 

PROCEDURE  FOR  SOLVING  PROBLEM  A 

(a)  Let  fp  be  the  flow  pattern  that  is  defined  by  f(u)  =  0  for  all  u  e  U. 
This  is  an  ideal  x^  -•x^  flow  pattern  of  value  zero. 

(b)  In  general,  having  determined  an  ideal  flow  pattern  f.  of  value  i, 
apply  Algorithm  1  and  find  a  simple  x^^  -♦x^  path  P.  such  that  e(P.;fp  is 
minimum. 

(c)  If  e(Pj;fp  =  «>  then  i  =  M,  i.e. ,  there  are  no  feasible  x^^  -*x^ 
flow  patterns  having  larger  values.  So  (fp,  f^^,  .  .  . ,  L  }  is  a  complete 
solution  to  Problem  A. 

(c’)  If  e(P.;f)  <  «  then  =  L  +  g.  is  an  ideal  x^^  -*x^  flow  pattern  of 
value  i  +  1,  where  g.  is  the  elementary  path  flow  corresponding  to  P^. 

One  can  now  replace  L  with  and  repeat  b. 

In  order  to  accelerate  the  process  of  finding  the  solution  set 
{  fp,  fj^,  .  .  it  is  desirable  to  increase  the  value  of  successive 
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patterns  by  more  than  one  unit  if  possible.  The  following  result  provides 
a  mechanism  for  doing  so. 


Lemma  7 

Let  f.  be  an  ideal,  nonmaximal  pattern  in  (X,  U)  having 

value  i,  and  let  be  a  simple  x^  -►x^  path  in  (X,  W)  that  minimizes 
e(P;f. ).  For  each  arc  we  P.  define  Q  as  follows: 

If  weU, 


Q  _  c(w)  -  f(w)  if  f(w)  >  0 
^  ' 


if  f(w)  <  0 


If  w  =  u'eU’, 


f  f(u)  -  b(u)  if  f(u)  <  0 
\f(u)  iff(u)>0 


Finally,  let  Q  =  p  Then  for  1  <  k  <  Q, 

f. . ,  =  f.  +  kg. 
i+k  1 

is  an  ideal  x^^  flow  pattern  of  value  i+k,  where  g^  is  the  elementary 
path  flow  corresponding  to  P.  and  (kg.  )(u)  =  k  •  g;(u)  for  all  ueU. 

Proof,  hi  order  to  prove  this  lemma  the  following  subsidiary  result 
will  first  be  established:  If  f.,  and  f.^g  are  ideal  flow  patterns 
having  values  i,  i+1,  and  i+2  then 

According  to  Lemma  1,  can  be  expressed  as 

where  g^  and  are  elementary  path  flows,  h  is  an  appropriate  zero 
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valued  flow  ( corresponding  to  a  sum  of  cycle  flows )  and  g- ,  g„ ,  and  h  are 

1  £i 

conformal.  Moreover  from  Lemma  2  it  follows  that  f.  +  g,  +  g„  and  f.  +  h 
are  feasible.  From  Lemma  3  one  obtains  the  relation 

T(f.^2)  -  T(f.  +  gj  +  gg)  >  T(f.  +  h)  -  T(f. )  . 

Since  L  is  ideal  the  right-hand  side  is  nonnegative.  Hence 

T((i^.2)>T(fj+gj+g2). 

But  since  f.,  „  is  ideal,  and  f.  +  g,  +  g  is  also  feasible  and  has  value  i+  2, 
equality  holds  and  f.  +  g,  +  gn  is  also  ideal.  Applying  Lemma  3  again 

11a 

yields  the  relation 

T(fi  +  g^  +  g2)  -  T(f.  +  g^)  >  T(f.  +  g2)  -  T(f. )  . 

Since  f is  ideal,  T(  f .  +  g^^ )  >  T(  f^^ )  and  T(  f .  +  g2 )  >  T(  f ) .  Hence 
T(fi+  g^  +  g^)  -  T(fj^^)  >T(f.^^)  -  T(f.)  . 

But  the  leftmost  term  is  equal  to  T(f.^2),  and  hence 

Now  let  f.  be  an  ideal  flow  pattern  of  value  i,  and  again  assume  that 
patterns  with  values  i+1  and  i  +  2  exist;  i.  e. ,  assume  that  M  >  i+  2.  Let 
R  be  a  simple  -♦x^  path  minimizing  e(P;f. ),  so  that  L  +  g^^  is  ideal, 
where  g_  is  the  corresponding  elementary  path  flow..  Also,  let  S  be  a 
simple  path  that  minimizes  e(P;f.  gj^)»  so  that  L  +  gj^  +  gg  is  ideal. 
From  the  inequality  established  earlier  in  the  proof  e(S;f.  +  gj^)  >  e(  R;f. ). 
Now  if  R  had  been  such  that  e(  R;f.  +  gj^)  =  e(  R;L  )  then  L  +  2gj^  would 
necessarily  also  be  ideal.  Similarly  if  e(R;f.  +  2gj^)  =  e(R;L  )  then 
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+  3gj^  will  be  ideal.  In  general,  1  +  kg^^  will  be  ideal  so  long  as 
e[R;f.  +  (k-  l)gj^]  =e(R;f.)  . 

Since  the  effective  cost  of  a  path  is  the  sum  of  the  effective  costs  of 
its  component  arcs,  f.  +  kg„  will  be  ideal  as  long  as 

1  XV 

e[w;f.  +  (k  ~  l)gj^]  =  e(w;f.) 

for  all  wcR.  If  w€  U  and  f(w)  >  0  each  successive  unit  of  flow  will  cost 
aj^(w)  until  w  is  saturated,  i.e.,  until  flow  equals. c(w).  Then  the  effec¬ 
tive  cost  changes  to  «>.  So  the  effective  cost  for  arc  w  changes  when 
k  =  c(w)  -  f(w).  If  w£  U  and  f(w)  <  0  each  unit  of  flow  added  in  the 
direction  associated  with  w  will  cancel  a  unit  of  flow,  resulting  in  a 
reduction  in  cost  of  -a2(w).  This  will  continue  until  all  flow  is  canceled, 
at  which  time  the  effective  cost  of  w  changes  from  -a2(w)  to  aj^(w).  So 
in  this  case  the  effective  cost  changes  when  k  =  -f(w).  The  two  cases 
that  arise  when  wc  U'  are  very  similar,  and  the  detailed  reasoning  is 
omitted.  In  all  cases  the  value  of  Q  given  in  the  statement  of  the  lemma 
is  the  smallest  value  of  k  such  that  e(w;f  +  kgj^)  f  e(w;L).  So  if  Q  = 
w€^R  ^w  ^  largest  integer  k  such  that 

Tlf.  +  kgj^]  -  Tlf.  +  (k-  l)gj^]  =Tlf.+  gj^]  -  Tlf.]  . 

Hence  f.  +  kg.^  is  ideal,  for  k  =  1,  2,  .  .  . ,  Q.  This  completes  the  proof. 

It  follows  from  the  fact  that  e(P^;f )  <  <®  that  k  >  1.  In  many 
applications  k  will  be  large  if  the  capacity  is  large,  and  the 

number  of  repetitions  of  step  b  of  the  above  procedure  will  be  greatly 
reduced. 
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An  algorithm  for  finding  that  minimizes  e(P;f. )  is  obtained  by 
treating  e(w;f. )  as  the  distance  from  x  to  y,  where  we  W  and  w  =  (x,y). 
The  problem  of  finding  P.  is  then  one  of  finding  the  shortest  distance 
from  to  through  the  network  (X,  W).  Efficient  algorithms  for 
finding  shortest  paths  exist,  provided  the  distances  are  such  that  no 
cycle  has  negative  total  length.  Fortunately  this  is  the  case.  The 
following  result  can  be  established; 

Lemma  8 

A  necessary  and  sufficient  condition  that  f  be  an  ideal  x^^  -*x^  flow 

pattern  in  (X,  U)  is  that  e(C;f)  >  0  for  every  simple  cycle  C  in  (X,  W). 

Proof.  If  f  is  ideal  then  certainly  e(  C;f )  >  0  for  all  simple  cycles  C. 

For  if  e(C;f)  <  0  for  some  C  then  from  Lemma  6  it  would  follow  that 

T(f  +  g^)  <  T(f),  contradicting  the  fact  thit  f  is  ideal. 

Now  suppose  that  e(  C;f )  >  0  for  all  simple  cycles  C.  Let  g  be  any 

Xj^  flow  pattern  having  the  same  value  as  f.  Using  Lemma  1,  g  can 

be  written  as  f  +  h.  +  .  .  .  +  h  ,  where  the  h.'s  are  conformal  cycle 
1  n  1 

flows  corresponding  to  appropriate  cycles  C, ,  C„,  .  .  .,  C  .  Repeated 
application  of  Lemma  3  yields 

T(f+  h^  +  .  .  .  +  h^)  -  T(f+  h^  +  .  .  .  +  h^_j)  >T(f+  h^)  -  T(f) 

T(f+h^+  .  .  .  +  h^_^)  -  T(f+ h^+  .  .  .  + Vl^  ■ 

• 

TCf+h^  +  hg)  -  T(f+hj)  >T(f+  hg)  -  T(f) 

T(f+  h^)  -  T(f)  >0  . 
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But  the  right-hand  side  of  each  inequality  is  >  0  by  assumption. 

Adding  the  left  sides  yields 

T(f  +  h^  +  .  .  .  +  h^)  -  T(f)  >  0  . 

Since  f  +  h^^  +  .  .  .  +  h^  =  g  it  follows  that  T(g)  >  T(f).  Since  g 
was  arbitrary  this  implies  that  f  is  ideal.  The  proof  is  complete. 

In  order  to  describe  a  specific  algorithm  having  the  characteristics 
of  Algorithm  1  suppose  that  the  arcs  of  W  are  arranged  in  a  fixed  sequence 
w. ,  w_,  .  .  .,  w  .  (It  is  particularly  useful  to  let  u  and  u' be  adjacent 
terms,  since  one  uses  the  same  basic  data  when  considering  u  and  u'.  ) 

For  convenience  let  e^  denote  e(Wj;f)  where  f  is  the  ideal  flow  pattern 
under  consideration.  The  following  procedure  has  the  desired 
characteristics. 

Algorithm  1 

(a)  Associate  numbers,  termed  labels,  with  vertices,  denoted  by 
V(x.)  as  follows: 

V(x^)  =  0 

V(xp  =  "0  for  i  =  2,  3,  .  .  .,  n  . 

(b)  Starting  with  w^^  take  each  arc  in  turn  and  do  the  following: 

(i)  If  w^  =  (Xj.Xj^)  and  V(Xj^)  >  V(Xj )  +  e.  replace  V(Xj^)  by  the 
smaller  quantity  V(x. )  +  e^  and  record  w.  as  the  "  approach  arc  "  asso¬ 
ciated  with  vertex  Xj^.  (If  this  is  not  the  first  iteration,  w.  supersedes 
any  approach  arc  previously  associated  with  x,^. ) 

(ii)  IfV(x^) 
arc  next. 


<  V(Xj)  +  e.  do  not  modify  V(Xj^),  but  consider 
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When  i  =  m  proceed  to  step  c. 

(c)  If,  in  the  course  of  considering  w, ,  w„ . w  in  step  b,  at 

X  z  m 

least  one  label  was  modified,  repeat  step  b. 

(c')  If  no  vertex  was  relabeled,  stop.  •  At  this  time  V(x^)  equals 

min  e(P;f)  taken  over  all  simple  x^  -•x^  paths  P.  If  V(x^)  =  «,  f  is  a 

maximum-valued  flow.  If  V(x^)  <  «>  trace  backward  along  the  approach 

arc  associated  with  x  to  its  initial  vertex  x  .  Trace  backward  along 

n  p 

X  's  approach  arc  to  x  and  continue  until  x.  is  reached.  This  will 
P  ^ 

yield  a  simple  x^^  -•x^  path  (traced  backward),  minimizing  e(P;f ). 

The  remainder  of  this  section  provides  a  proof  of  the  fact  that  this 
algorithm  will  always  produce  the  desired  information  in  a  finite  number 
of  steps  regardless  of  the  network  (X,  U)  and  the  ideal  flow  f  under  con¬ 
sideration.  Some  of  the  properties  of  this  algorithm  will  be  needed  in  a 
more  general  form  in  connection  with  Problem  C.  For  this  reason  the 
assumption  that  f  is  ideal  is  not  made  in  establishing  a  number  of 
intermediate  results. 

One  "iteration"  of  the  algorithm  will  mean  one  complete  pass  through 

the  entire  sets  of  arcs  {  w, ,  w„,  .  .  . ,  w^  },  relabeling  a  vertex  whenever 

X  z  in 

appropriate.  From  the  time  that  a  vertex  other  than  x^^  first  attains  a 
finite  label,  it  has  associated  with  it  an  approach  arc  leading  from 
another  vertex  that  must  also  have  a  finite  label.  The  approach  arc  asso¬ 
ciated  with  a  vertex  may  change  from  time  to  time.  Since  the  approach 
arc  involved  in  assigning  a  new  label  supersedes  the  previous  one  (if  any) 
identified  with  the  vertex  being  relabeled,  at  no  time  does  a  vertex  have 
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more  than  one  approach  arc.  Summarizing,  if  V(x)  =  «,  x  has  no  approach 
arc.  If  V(x)  <  «,  X  has  precisely  one  approach  arc,  but  it  is  not  necessarily 
the  same  arc  throughout  the  procedure. 

Suppose  that  the  algorithm  is  started  and  then  interrupted  at  some  time, 
and  that  V(x)  <  «»  at  this  time,  where  x  f  x^^.  Then  x  must  have  associated 
with  it  an  approach  arc,  say  =  (yj^.x).  Similarly  V(yj^)  <  «>,  so  y^^  has 
an  approach  arc  =  (y2»yj^).  Repeating  this  reasoning  Vg  s  (yg.yg)  is 
obtained,  etc.  In  this  process  of  tracing  back  along  approach  arcs  one  of 
three  situations  may  arise.  These  are: 

(a)  Vertex  Xj^  may  be  reached  via  approach  arcs  v^,  Vg,  .  .  .,  v^ 

(i.e. ,  Xj^  is  the  initial  vertex  of  v^),  and  it  may  be  impossible  to  trace 

back  along  another  approach  arc  because  Xj^  has  no  associated  approach 

arc.  This  will  be  the  case  if  V(x)  =  0,  so  that  the  original  label  of  x^^ 

has  not  been  altered.  In  this  case  the  procedure  has  traced  (backward)  a 

simple  X. -•x  path  of  approach  arcs  (v  ,  v  v-,  v,  }.  (This 

X  n  m  in*x  ^  i 

path  is  necessarily  simple  because  no  vertex  has  two  approach  arcs 
leading  to  it. ) 

(b)  The  tracing  procedure  may  lead  to  Xj^,  but  x^^  may  have  been 
relabeled.  In  this  case  one  can  continue  to  trace  backward,  and  must 
inevitably  return  to  some  vertex  reached  earlier  since  each  step  leads 

to  another  vertex  with  a  finite  label  and  every  such  vertex  has  an  approach 
arc  leading  to  still  another  vertex.  Once  a  vertex  is  repeated  in  this 
tracing  process,  a  cycle  of  approach  arcs  has  been  traversed  (backward) 
and  further  tracing  will  merely  duplicate  the  route  already  covered.  In 
this  case  the  set  of  approach  arcs  discovered  by  tracing  takes  the  form  of 
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a  cycle  together  with  a  path  from  the  repeated  vertex  to  x.  ( If  the  vertex 
repeated  is  x  itself  the  set  of  arcs  is  simply  a  cycle.  ) 

(c)  The  tracing  procedure  may  repeat  a  vertex  without  ever  reaching 
Xj^.  In  this  case  the  form  of  the  set  of  approach  arcs  traced  is  the  same 
as  in  situation  b,  being  either  a  cycle  together  with  a  path  from  some 
vertex  of  the  cycle  to  x,  or  else  simply  a  cycle  containing  x.  The  only 
difference  is  that  x^^  is  not  incident  with  any  of  the  arcs  traced. 

In  the  remainder  of  this  section  for  x  f  Xj^,  A(x)  will  denote  the  set 
of  approach  arcs  traced  backward  from  x,  assuming  that  V(x)  <  i.e. , 
that  X  has  an  approach  arc.  The  set  A(x)  is  called  the  approach  set 
associated  with  x.  It  will  be  shown  that  if  the  flow  pattern  f  under  con¬ 
sideration  is  ideal,  case  a  always  holds,  so  that  A(x)  is  always  a  simple 
-•x  path.  If  f  is  not  ideal,  sooner  or  later  case  b  or  case  c  arises  for 
at  least  one  vertex  x.  The  important  fact  then  is  that  there  is  a  cycle  of 
approach  arcs  and  that  such  a  cycle  can  be  used  to  reduce  the  cost  of  f. 

The  symbol  A  will  denote  the  total  set  of  approach  arcs  associated 
with  vertices  at  a  given  time  during  execution  of  the  algorithm. 

Lemma  9 

If  Algorithm  1  is  applied  to  a  feasible  x^^  -►x^  flow  pattern  and  if  the 
procedure  is  halted  at  any  time  then,  if  wcA  at  this  time  and  if  w  s  (x,  y), 

V(y)  >  V(x)+  e(w;f) 

Proof.  Note  first  that  in  this  and  in  subsequent  lemmas  that 
investigate  the  status  at  an  arbitrary  time  during  the  algorithm,  the  labels 
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and  approach  arcs  referred  to  are  those  in  effect  at  the  time  of  interruption 
unless  the  contrary  is  specified. 

Since  w  is  the  approach  arc  associated  with  y  at  this  time  it  follows 
that  V(y)  =  V*(x)  +  e(w;f)  where  V*(x)  is  the  label  that  was  associated 
with  X  at  the  time  y  attained  its  present  label.  But  V*(x)  >  V(x)  since 
the  procedure  never  increases  labels.  It  follows  that  V(y)  >  V(x)  +  e(w;f), 
which  completes  the  proof. 

Lemma  10 

If  Algorithm  1  is  applied  to  a  feasible  -►x^  flow  pattern  and  if  P  is 
any  simple  x^^  -•y  path  having  k  arcs  then  after  k  or  fewer  iterations  of 
the  algorithm  the  following  relation  holds: 

V(y><e(P;f)  . 

Proof.  LetP  =  {wj^,  Wg . Wj^}  where  w^  s  (x^,y^ ),  W2S(y^,y2), 

.  .  .,  Wr  s!  (yR_^,,y).  Let  V^(2)  denote  the  label  of  vertex  z  after  t  iterations, 
with  the  convention  that  Vq(z  )  denotes  the  original  label  of  z.  Then  clearly 

<  Vo(Xi)  +  e(Wj;f) 

since  every  arc,  including  w^,  is  considered  as  a  candidate  for  relabeling 
Yy  Similarly, 

Vl)''  • 

By  repeating  this  reasoning  this  relation  is  ultimately  obtained: 

VR(y)  <  Vq(Xj)  +  sj  e(w.;f)  =  e(P;f)  . 

The  second  relation  holds  since  Vq(x^  )  =  0.  If  the  algorithm  terminates 
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after  fewer  than  k  iterations  then  at  the  time  of  its  termination 
V(yj^)  <  V(x^)  +  e(w^;f),  .  .  V(yj^)  <  V(yj^_^)  +  e(Wj^;f).  Since 
V(Xj^)  <  0  this  again  implies  that  V(yj^)  <  e(P;f ).  This  completes  the 
proof. 

Lemma  11 

If  Algorithm  1  is  applied  to  a  feasible  x^  -*x^  flow  pattern  f  and  if  at 
some  time  the  set  of  approach  arcs  A  contains  a  cycle  then  f  is  not  ideal. 

Proof.  Suppose  that  C  =  {  Wj,  w^,  .  .  . ,  Wj^}  is  a  set  of  approach 
arcs  that  forms  a  cycle.  (Note  that  C  must  be  a  simple  cycle  because  no 
two  arcs  of  A  ever  terminate  at  the  same  vertex. )  Assume  that  Wj^  s  (y^^,  yg ), 
•  •  • »  Wr  —  (yj^iyj^)*  By  applying  Lemma  9  the  following  inequalities  are 
obtained: 

VCyg)  >V(y^)+  e(w^;f) 

V(y3)  >  V(y2)  +  e(W2;f) 

• 

« 

V(y^)  >  V(yj^)+  e(wj^;f)  . 
k 

It  follows  that  0  >  e(w^;f )  =  e(C;f ).  Now  let  y^  be  the  first  vertex  of 
the  set  {  y^^,  y2,  .  .  . ,  y^^}  which  attained  its  present  label.  Then 
V(yj)  =  V*(yj  +  e(w^;f),  where  V*(y^_^)  was  the  label  of  y^  at  that 
time.  By  the  way  in  which  y^  was  chosen  Vfy^  <  V*(yj_j^).  It  follows 
that  one  of  the  above  inequalities  is  strict,  so  that  one  obtains  the  sharper 
relation  0  >  e(C;f).  But  then  it  follows  from  Lemma  8  that  f  cannot  be 
ideal.  This  completes  the  proof. 
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Lemma  12 


If  Algorithm  1  is  applied  to  a  feasible  -*x^  flow  pattern  f  and  if 
V(x^ )  <  0  at  some  stage  then  from  this  time  on  A  contains  a  simple  cycle. 

Proof.  Assume  the  algorithm  is  stopped  some  time  after  V(Xj^) 
becomes  negative,  with  A  denoting  the  set  of  approach  arcs  at  the  time  of 
stopping.  Let  w^  be  the  approach  arc  associated  with  x^^,  and  let  be 
the  vertex  such  that  w^^  =  (y^^.x^).  Similarly  let  Wg  be  the  approach  arc 
associated  with  y^,  and  let  y^  be  such  that  W2  =  (y2,  y^  )•  can  con¬ 
tinue  to  trace  backward  obtaining  (wg,  w^,  .  .  .  }  and  {yg,  y^,  .  .  .  }. 

Let  k  be  the  smallest  integer  such  that  Wj^  s  (yj^.yj^  )i  and  such  that  yj^ 
coincides  with  or  with  y^  for  some  j  <  k.  (A  repetition  must  occur 
eventually,  since  there  are  only  a  finite  number  of  vertices.  )  If  yj^  =  x^^ 

then  {w^,  Wg . Wj^}  is  a  simple  cycle.  If  yj^  =  yj  then  w-^g* 

.  .  . ,  Wj^}  is  a  simple  cycle.  In  any  case  A  contains  a  simple  cycle. 

This  completes  the  proof. 

The  following  result  is  an  immediate  consequence  of  the  last  two 
lemmas: 

Lemma  13 

If  Algorithm  1  is  applied  to  a  feasible  flow  pattern  f  and  if  at 

some  time  V(Xj^)  <  0  then  f  is  not  ideal. 

Lemma  10  established  the  fact  that  V(y)  is  less  than  or  equal  to  the 
total  effective  cost  e  (P;f)  of  any  x^  -•y  path  P  after  a  sufficient  number 
of  iterations.  Lemma  14  establishes  a  reverse  inequality.  The  two 
results  together  are  then  used  to  show  that  after  a  sufficient  number  of 
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iterations  e(P;f)  is  precisely  equal  to  V(y),  provided  that  f  is  ideal  and 
that  y  ever  attains  a  finite  label. 

Lemma  14 

If  Algorithm  1  is  applied  to  a  feasible  flow  pattern  f  and  if  the 

process  is  interrupted  at  any  time  then,  if  A(y)  is  a  simple  x^^  -•y  path 
P  for  some  vertex  y  =  x^^,  the  inequality  V(y)  >  e(P;f)  must  hold. 

Proof.  Suppose  P  =  {wj^,  w^,  .  .  .,  Wj^}  where  s  (x^.y^^), 

Wg  ®  (yi.yg).  .  .  ..  S  (yj^_j^,y).  it  follows  from  Lemma  9  that 

V(yi)  >V(x^)+  e(w^;f) 

V(y2)  >  V(yj^)+  e(W2;f) 


V(y)  >  (\_i)  +  e(Wj^;f)  . 

Summing,  and  using  the  fact  that  V(x^)  =  0  [because  A(y)  is  a  simple 
path]  V(y)  >  e(P;f)  is  obtained.  This  completes  <'he  proof. 

Theorem  15 

If  Algorithm  1  is  applied  to  a  nonmaximal  ideal  x^^  -•x^  flow  pattern 
in  a  Tjq)e  A  network  having  n  vertices  the  algorithm  will  terminate  ( cease 
to  reduce  any  labels)  after  n  or  fewer  iterations  or  complete  "passes" 
through  the  set  of  arcs.  On  termination,  if  V(y)  <  *,  A(y)  will  be  a 
simple  -*y  path  that  minimizes  e{P;f )  over  all  such  paths.  In 
particular  A(x^)  will  be  an  optimal  path. 
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Proof.  Since  f  is  ideal.  Lemma  11  states  that  at  no  time  will  the  set 


of  approach  arcs  contain  a  cycle.  So  for  every  y  ^  Xj^,  either  V(y)  =  «> 
or  else  A(y)  is  a  simple  -*y  path  P.  In  the  latter  case  V(y)  >  e(P;f) 
according  to  Lemma  14.  So  V(y)  is  always  as  large  as  the  effective  cost 
of  some  simple  x^  -*y  path.  But  if  P^  is  an  x^  ^y  path  that  minimizes 
e(P;f ),  and  if  P^  has  k  arcs  then  V(y)  <  e(Py.;f)  after  k  or  fewer  iterations, 
according  to  Lemma  10.  Strict  inequality  cannot  hold,  so  V(y)  =  e(Py;f) 
after  k  or  fewer  iterations.  But  k  <  n  because  a  simple  path  cannot  have 
more  than  n  arcs.  So  those  vertices  other  than  x^  that  ever  have  finite 
labels  achieve  their  smallest  labels  after  n  or  fewer  iterations.  As  for 
Xj^  itself  Lemma  13  asserts  that  it  is  never  relabeled.  It  remains  to  show 
that  if  V(y)  <  <»  on  termination  then  A(y)  is  a  simple  x^  -►y  path  of  mini¬ 
mum  effective  cost.  [  Note  that  P^  was  any  effectively  "  cheapest "  path 
in  the  above  argument.  A(y)  may  be  a  different  one.  ]  Let  P^^  denote 
A(y)  on  termination.  Since  V(y)  =  e(Py;f),  and  since  e(Pj^;f)  >  e(Py;f) 
because  P^  is  optimal  it  follows  from  Lemma  10  that  e(Pj^;f)  =  e(Py;f), 
so  Pj^  also  minimizes  e(P;f).  In  particular  if  y  =  x^,  V(x^)  <  ■»  on 
termination  because  f  is  not  maximal,  by  assumption.  So  the  above 
remarks  apply,  and  A(x^)  will  be  a  simple  x^^  -•x^  path  minimizing 
e(P;f ).  This  completes  the  proof. 

If  Algorithm  1  is  applied  to  a  maximal  ideal  x^^  ^x^  flow  pattern  the 
only  difference  is  that,  on  termination  (which  still  must  occur  after  n  or 
fewer  iterations),  V(x^)  =  «.  Hence  if  the  algorithm  is  applied  to  any 
ideal  pattern  the  procedure  will  terminate  after  n  or  fewer  iterations.  If 
V(x^)  =  «,  f  is  maximal.  If  V(x^)  <  «  then  A(x^)  is  an  x^^  -•x^  path 
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minimizing  e(P;f).  Hence  Algorithm  1  has  the  properties  required  for 
the  procedure  given  earlier  to  solve  Problem  A.  (Several  of  the  lemmas 
used  to  establish  this  fact  will  also  be  used  to  devise  a  procedure  for 
solving  Problem  C. ) 
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SOLUTION  OF  PROBLEM  B 


This  section  describes  a  complete  solution  to  Problem  B.  Recall 
that  a  Type  B  network  is  a  weighted  and  capacitated  network  (X,U)  such 
that  b(u)  >  0  for  at  least  one  arc  uc  U.  Moreover  two  vertices  and 
are  desig^iated  as  the  source  and  sink  of  all  flow  patterns  under  consider¬ 
ation.  The  problem  is  that  of  finding  an  ideal  x^^  "•Xn  pattern  for 
each  value  realized  by  feasible  Xj^  -►x^  patterns — with  the  possibility  that 
there  may  be  no  such  values  in  the  given  network. 

In  essence  the  procedure  given  here  for  solving  Problem  B  consists 
of:  (a)  associating  with  the  Type  B  network  a  related  Type  A  network 
(X,  U),  (b)  solving  Problem  A  for  the  related  network,  and  (c)  interpreting 
the  resulting  ideal  flow  patterns  as  the  solution  to  the  Type  B  problem. 

The  construction  of  the  associated  network  (X,  U)  employed  here  is  a 

2 

modification  of  a  construction  used  by  Berge,  and  the  procedure  for 
finding  a  complete  set  of  ideal  Type  B  flows  is  an  extension  of  a  procedure^ 
indicated  in  that  reference,  for  determining  a  feasible  Type  B  flow  pattern 
of  least  cost. 
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Given  a  Type  B  network  (X,  U);  construct  a  Type  A  network  (X,  U)  as 
follows: 

(a)  Augment  X  by  adding  two  vertices  and  x^.  ( This  notation  is 
adopted  because  these  vertices  will  constitute  the  source  and  sink  of  flow 
patterns  in  the  augmented  network. )  Let  X  denote  X  together  with  x^^  and 

Xg. 

(b)  Let  S  denote  the  set  of  arcs  of  U  such  that  b(u)  >0.  (The  set  S 
occurs  frequently  in  this  section  and  always  has  the  above  meaning.  )  If 
ucS  and  u  s  fx,y),  construct  two  additional  arcs  u*  and  u**  where 

u*  s  (Xj^.y)  and  u**  =  (x,x^).  Let  U“  and  U**  denote  the  set  of  arcs  of 
type  u*  and  u**  respectively.  Each  of  these  sets  is  in  one-to-one 
correspondence  with  S. 

(c)  Construct  two  additional  arcs  w*  and  w“*,  where  w*  =  (x^.x^^) 
and  w**  S'  (x^,x^).  Let  U  denote  the  totality  of  arcs  of  U,  U*,  and  U**, 
plus  the  two  arcs  w*  and  w**.  (Figure  5  is  an  example  intended  to  clarify 
the  above  construction.  In  this  example  the  heavy-line  arcs  constitute 
the  set  S. ) 

(d)  .  Lower  and  upper  arc  bounds  denoted  by  b(u)  and  c(u)  for  ucU 
are  defined  as  follows: 

b(u)  =  b(u)  and  c(u)  =  c(u)  if  ucU-S.* 
b(u)  =  0  and  c(u)  =  c(u)  -  b(u)  if  ucS. 
b(u* )  =  0  and  c(u*)  =  b(u)  for  u*  c  U*. 

’'‘U-S  denotes  the  arcs  of.  U  that  are  not  in  S,  i.e. ,  those  ue  U  such 
thatb(u)<0. 
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o.  Givan  nafwerk  (X,  U)  t 


b.  Assoc  iatad  natwork  (X,  U)'|’ 


Fig.  5 — Ralotiofl  batwaan  Natworks  (X,  U)  and  (X,  U) 
tBold  area  indicota  lha  sat  S. 


b(u")  =  0  and  c(u**)  =  b(u)  for  u**eU". 
b(w* )  =  b(w**  )  =  0 
c(w*)  =  c(w**)  =  K^. 

Here  is  any  integer  known  to  exceed  the  value  of  every  feasible  -•x^ 
flow  pattern  in  (X,U). 

(e)  Unit  costs  denoted  by  aj^(u)  and  defined  as  follows: 

a^(u)  »  aj^(u)  and  a2(u)  =  a2(u)  for  ucU. 
aj^(u* )  =  a2(u* )  =  0  for  u*  e  U*. 
aj^(u“  )  =  a2(u** )  =  0  for  u“  e  U**. 
aj^(W)  =  a^(w**)  =  K2. 
a2(W)  =a2(W)  =  0. 

Here  K„  is  an  integer  that  exceeds  S..max{a,  (u),  a  (u)}.  (The  reasons 
for  specifying  and  K2  in  this  manner  will  be  clarified  later.  ) 

If  Xj^  and  x^  are  designated  as  source  and  sink  in  the  augmented 
network  then  (X,  U)  together  with  functions  b,  c,  a^^,  and  define  a  Type 
A  network.  Whenever  feasible  or  ideal  flow  patterns  in  (X,  U)  are  men¬ 
tioned  these  are  intended  to  be  x^  patterns  that  are  feasible  or  ideal 
relative  to  b,  c,  a^^,  and  a2.  Patterns  in  (X,  U)  are  always  related  to  b, 
c,  a^  and  a2.  Flow  patterns  in  (X,  U)  will  be  denoted  by  symbols  such 
as  7 and  g  to  further  distinguish  them  from  flows  in  (X,  U). 

Let  7  be  a  feasible  flow  pattern  in  (X,  U)  which  saturates  U*  and  U**, 
i.e. ,  one  such  that  f(u* )  =  c(u*  )  for  all  u*  c  U*  and  f(u** )  =  c(u”  ) 
for  all  u”  €  U**.  Such  patterns,  and  only  such,  are  called  transformable 
flow  patterns  in  (X,  U).  It  will  be  shown  that  transformable  flow  patterns 
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in  (X,  U)  correspond  to  feasible  patterns  in  (X,  U)  and  that  ideal 
transformable  patterns  in  (X,  U)  correspond  to  ideal  patterns  in  (X,  U). 

Let  f  be  an  arbitrary  transformable  flow  pattern  in  (X,  U).  A 
transformation  <f)  will  now  be  introduced  that  maps  f  into  a  flow  pattern  f 
in  (X,  U).  Specifically,  for  every  uc  U  define  f(u)  as  follows: 


f(u) 


f(u)  ifucU-S 

f(u)  +  b(u)  if  ucS 


The  pattern  f  related  to  f  in  this  way  will  be  denoted  sometimes  by  <pi. 

Note  that  if  f^^  and  are  different  transformable  flow  patterns  in  (X,  U) 
then  0fj^  and  will  be  different  functions  on  U.  For  if  f^^  and  coincide 
on  U  they  will  coincide  also  on  U*,  U**,  w*,  and  w**.  (All  transformable 
flow  patterns  have  the  same  values  on  arcs  of  U“  and  U** ,  so  that  it  is 
only  necessary  to  consider  w*  and  w** .  But  f^^  and  fg  must  be  such  that 
input  equals  output  at  and  Since  the  values  of  arc  flow  are  the  same 
for  all  arcs  incident  with  Xj^,  except  possibly  w*,  they  must  agree  also  on 
w* .  A  similar  consideration  applies  for  w** .  )  Thus  if  f  and  are 
distinct  they  must  differ  on  an  arc  of  U,  so  that  f^^  =  will  not  coincide 
with  fg  =  Expressed  differently,  <f>  is  a  one-to-one  mapping  of  the 
set  of  transformable  x^  -•x^  flow  patterns  in  (X,  U)  onto  a  certain  set  of 
flow  patterns  in  (X,  U).  This  latter  set  is  in  fact  the  set  of  all  feasible 
Xj^  -•x^  flow  patterns  in  (X,  U),  as  will  be  seen  presently. 


Lemma  16 

If  f  is  a  transformable  Xj^  -^x^  flow  pattern  of  value  k  in  (X,  U)  then 
f  =  <^is  a  feasible  x^  -•x^  flow  pattern  of  value  k  -  Egb(u)  in  (X,  U). 
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Proof.  Note  first  that  each  arc  flow  is  feasible.  For  if  ue  U-S, 
f(u)  =  f(u),  b(u)  =  b(u)  and  c(u)  =  c(u),  so  that  the  feasibility  of  f(u) 
follows  from  that  of  f(u).  If  ucS  then  f(u)  =  f(u)  +  b(u),  b(u)  =  0,  and 
c(u)  =  c(u)  -  b(u).  Since  f  is  feasible  0<  f(u)  <  c(u).  Consequently 
.5  .5  c(u),  so  f(u)  is  feasible. 

It  will  be  shown  next  that  if  x«X  and  x  is  neither  x,  nor  x  the  net 

1  n 

output  at  X  is  zero.  The  summation  2u(-*x)  written  as 

2j^f(u)  +  S2f(u)  where  ranges  over  arcs  of  U  n  U('-»x)*  and  Sg 
ranges  over  U*n  U(-*x).  For  every  ucS  fl  U(-*x),  contains  a  term 
f(u),  and  Sg  contains  a  term  c(u*  )  =  b(u)  which  together  equal 
f(u)  +  b(u)  =  f(u).  For  every  uc(U-S)  n  U("*'X).  contains  a  term 
f(u),  and  has  no  corresponding  term,  since  b(u)  <  0  for  arcs  of  U-S 

A 

and  no  corresponding  arc  u*  exists.  Also,  f(u)  =  f(u)  in  this  case.  It 

follows  that  S.,,  ,  f(u)  =  2)7:. _  V  f(u)  since  each  term  of  the  first  sum- 

U('^s)  '  U("*x) 

mation  can  be  paired  with  one  or  two  terms  of  the  second  summation, 
depending  on  whether  uc  (U-S)  or  whether  ucS.  A  very  similar  argu¬ 
ment  shows  that  S...  =  2-=:.  Since  the  net  output  of  f 

U(X“*)  U\X~») 

at  X  is  zero  it  follows  that  the  net  output  of  f  at  x  is  zero.  Thus  f  is  either 
a  feasible  -•Xj^  flow  pattern  or  a  feasible  x^  -*Xj^  flow  pattern  in  (X,  U). 

It  will  next  be  shown  that  the  former  possibility  is  always  true.  (The 
latter  will  also  be  true  whenever  the  value  of  f  is  zero. ) 

Consider  now  vertex  x=Xj.  One  can  write  S-y^_^^jf(u)  =  Sj^f(u)  +  Z;2f(u)  + 
7(w* )  where  2  ^  and  Zg  have  the  same  meanings  as  before.  By  matching 

’^A  ri  B  denotes  the  set  of  elements  common  to  sets  A  and  B. 
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all  terms  of  and  with  terms  of  S^^_^^^f(u),  only  the  arc  flow  f(w* ) 
has  no  counterpart  in  f.  So  =  S— ^_^^^f(u)  -  f(w* ).  As 

before  S-  .  _*vf(u)=  S-^,  Since  net  output  of  f  at  x  is  zero 

it  follows  that  the  net  output  of  f  at  is  f  (w“ ).  So  f  is  a  feasible  x^^  -►x^ 
flow  pattern  of  value  f(w® )  in  (X, U).  Now  the  value  of  f  is  k  =  f(u)  + 
f(w‘’ )  and  Syo  f(u)  =  Sgb(u).  Hence  the  value  of  f  is  k  -  Sgb(u)  as 
asserted.  This  completes  the  proof. 

It  was  noted  above  that  no  flow  pattern  in  (X,  U)  is  the  transform  of 
more  than  one  transformable  flow  pattern  in  (X,  U).  On  the  other  hand 
it  can  readily  be  shown  that  an  arbitrary  feasible  x^  pattern  f  in 

(X,  U)  is  the  transform  of  some  transformable  flow  pattern  in  (X,  U).  To 
show  this,  define  a  transformation  ♦  as  follows:  ♦f  =  f,  where  f  is  the 
flow  pattern  in  (X,  U)  defined  by: 


f(u)  =  f(u)  if  ucU-S. 

f(u)  =  f(u)  -  b(u)  if  u€ S. 

f(u*)  =  b(u)  if  u*€U*. 

f(u“)  =b(u)  if  u*’cU”. 

f(w*  )  =  f(w**  )  =  k,  where  k  is  the  value  of  f. 

It  is  clear  that  f  is  a  feasible  x^^  flow  pattern  in  (X,  U).  ( The  fact 
that  net  output  is  zero  for  all  xcX  is  established  by  introducing  2^  and  Sg 
as  before  and  applying  essentially  the  same  argument  as  before. )  The 
value  of  7  is  k+  S_b(u)  where  k  is  the  value  of  f. 

The  transformation  ♦  is  the  inverse  of  (p,  satisfying  ^(p  f )  =  ^  =  f 
and  4){^)  =  <pJ=  f.  So  ♦will  be  denoted  by  <p~^.  Summarizing  the  above 
observations  gives  Theorem  17. 
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Theorem  17 


(pis  a  one-to-one  mapping  of  the  set  of  all  transformable  -*x^  flow 
patterns  in  (X,  U)  onto  the  set  of  all  feasible  flow  patterns  in  (X,  U). 

If  f  is  transformable  and  has  value  k  then  (pf  has  value  k  -  S  b(u).  The 

O 

transformation  ♦  defined  above  is  the  inverse  mapping  ^  of  all  feasible 

X,  -►X  flows  onto  all  transformable  x,  -*x  flows. 

In  in 

The  costs  as  well  as  the  values  of  a  transformable  x,  -*x  flow  pattern 

In 

f  and  its  image  f=  (pi  are  connected  by  a  simple  relation.  The  following 
result  establishes  this  relation. 

Theorem  18 

If  f  =  <pf  then 

T(f )  =  T(f )  -  2a^(w*  )f  (w* )  +  2ga^(u)b(u)  . 

Since  f(w* )  was  seen  earlier  to  be  the  value  of  f,  and  Sgaj^(u)b(u)  is  a 
constant,  this  means  that  T(f )  differs  from  T(f)  by  an  additive  constant 
and  a  term  proportional  to  the  value  of  f. 

Proof.  For  u€ U-S,  f(u)  =  f(u)  and  a.(u)  =  a.(u)  for  i  =  1  and  2,  so 
the  costs  of  f  and  f  are  the  same  on  U-S.  For  ue  S  the  cost  of  f  is 
Sgaj^(u)f(u)  and  that  of  Tis  Sgaj^(u)  lf(u)  -  b(u)],  so  that  the  cost  of  f 
exceeds  that  of  Tby  2gaj^(u)b(u)  on  S.  Since  f  is  only  defined  on  arcs  of 
S  and  U-S  the  total  cost  of  f  is  accounted  for  above.  Since  aj(u)  =  a2(u)  -  0 
for  ut  U*  or  u€  U**  the  only  additional  cost  of  f  is  associated  with  arc  flows 
7(w* )  and  T(w** ).  But  since  T{w*  )  =7(w**  )  >  0  for  transformable  flows, 
and  aj^(w* )  =  aj^(w** )  by  definition,  the  total  cost  in  w*  and  w"  is 
2aj^(w*  )7(w* ).  The  formula  given  in  the  statement  of  this  theorem  is 
simply  a  consolidation  of  the  above  relations. 
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As  a  corollary  to  the  preceding  theorem  it  may  be  noted  that  if  f  is  a 
feasible  -•x^  flow  pattern  of  value  k,  and  f  =  then 

T(f)  =  T(f)  -  2ki^J(w•  )  +  Sga^(u)b(u)  . 

so  that  the  costs  of  f  and  f  differ  by  a  constant  if  the  value  of  f  is  fixed. 
Therefore  to  find  a  minimum-cost  (i.e. ,  ideal)  Xj^-*x^  flow  pattern  of 
value  k  it  suffices  to  find  a  minimum-cost  transformable  flow 

pattern  f  such  that  the  value  of  ^f  is  k.  Since  the  value  of  f  exceeds  that 
of  f  =  (^f  by  £gb(u)  it  suffices  to  find  a  transformable  f  whose  value  is 
k  +  £gb(u)  and  such  that  no  other  transformable  flow  having  this  value 
has  smaller  cost. 

Nearly  all  the  information  required  for  a  general  solution  to  Problem  B 
is  now  at  hand.  The  remaining  gap  to  be  filled  is  to  show  that  every  non- 
transform  able  x^-*x^  flow  pattern  of  a  given  value  costs  more  than  a 
minimum-cost  transformable  x^-*x^flow  pattern,  provided  that  any  trans¬ 
formable  patterns  having  this  value  exist.  For  otherwise,  if  one  simply 
computes  ideal  Xj^-^x^  flows  for  all  feasible  values,  these  flows  may  not 
be  transformable  and  hence  will  not  yield  corresponding  ideal  Xj^-*x^ 
flows.  The  next  result  removes  this  danger,  hi  the  course  of  doing  so, 
the  reason  for  associating  a  large  unit  cost  (viz, ,  K^)  with  w*  and  w** 
will  be  clarified. 

Lemma  19 

If  any  transformable  flow  patterns  of  value  k  exist  then  every 

ideal  pattern  of  value  k  is  transformable. 
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Proof.  Assume  that  and  are  feasible  patterns  having 

the  same  value  k,  and  that  L  is  transformable  whereas  f_  is  not.  It  will 

^  a 

be  shown  that  f^  cannot  be  ideal.  Consider  the  zero-valued  flow  pattern 
fj^  -  fg.  By  Lemma  1  it  is  possible  to  write  “  ^2  ~  ^  ^  ‘  ^  ®n 

where  the  gJs  are  appropriate  conformal  elementary  cycle  flows.  Since 
fg  is  not  transformable  U*  and  U**  are  not  both  saturated  by  whereas 
f^  does  saturate  both  U°  and  U**.  Suppose  f^  fails  to  sa 'urate  U".  Then 
”  ^2  positive  flow  in  one  or  more  arcs  u'eU*,  zero  flow  in  the 

remaining  arcs  of  U* ,  and  negative  flow  (i.  e. ,  flow  toward  )  in  w* .  At 
least  one  of  the  elementary  cycle  flows,  say  g.,  must  be  such  that  g^(u® )  =  1 
for  some  u*  f  U®  and  g.(w®  )  =  -1,  since  the  sum  of  all  of  the  eyelet  ac¬ 
counts  for  all  flow  of  ~  ^2*  Lemma  2  it  follows  that  +  g^  is  a 

feasible  flow  pattern  whose  value  is  k.  It  will  now  be  shown  that 
T(f2  +  g. )  <  T(f2),  so  that  cannot  be  ideal.  Now  e(g.;f2)  is  the  sum  of 
terms  each  of  which  is  either  plus  or  minus  the  unit  cost  of  an  arc  of  the 
elementary  cycle  associated  with  g..  In  particular  one  of  the  terms  is 
-K2  where  was  the  unit  cost  of  w®,  i.e. ,  aj^(w®  ).  If  the  cycle  also 
contains  arc  w*®,  g^(w®® )  =  -1  for  the  same  reasoning  as  that  applied  to 
w®.  So  e(g.;f„)  =  -K  +  Q  or  -2K  +  Q  where  Q  is  the  sum  of  the  effective 

1  M  A  A 

lengths  of  arcs  of  U  involved  in  the  cycle.  ( The  cycle  also  contains  an  arc 
of  type  u*  and  possibly  one  of  type  u®®,  but  since  these  have  zero  unit  costs 
their  effective  lengths  are  zero  and  can  be  disregarded.  )  But  it  is  clear 
that  K2  >  Q  since  Kg  was  taken  to  be  an  integer  greater  than 
S^max  {aj^(u),  SgCu)}  and  hence  is  surely  greater  than  the  sum  of  the 
effective  lengths  relative  to  fg  of  selected  arcs  of  U.  It  follows  that 
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<  0  so  that  T(f2  +  gp  <  T(f2)  by  Lemma  6.  Hence  is  not  ideal. 

( If  U*  were  saturated  by  but  U**  were  not,  the  proof  would  be  essentially 
the  same. ) 

All  the  necessary  facts  that  lead  up  to  the  following  theorem,  whose 
validity  forms  the  basis  for  a  general  procedure  for  solving  Problem  B, 
have  now  been  established. 

Theorem  20 

Let  M  denote  the  capacity  of  (X,  U).  Either  there  are  no 

transformable  flow  patterns,  or  else  there  is  an  integer  m,  where 

0  <  m  <  M,  such  that  there  are  no  transformable  Xj^-*x^  flow  patterns  of 
value  p  for  0  <  p  <  m,  and  every  ideal  pattern  of  value  q,  where  m  <  q  <  M, 
is  transformable. 

Proof.  Suppose  there  exists  a  transformable  flow  pattern  of 

value  i,  where  i  <  M.  It  follows  from  Lemma  19  that  any  ideal  x,  -•x 
flow  pattern  of  value  i  is  necessarily  transformable.  Let  L  be  such  a 
pattern.  Since  f.  is  ideal  but  not  maximal  it  follows  from  Lemma  5  that 
an  ideal  flow  pattern  of  value  i+ 1  can  be  obtained  by  adding  an  appro¬ 
priate  elementary  path  flow  to  f..  But  since  every  arc  of  type  u*  and  u** 
is  already  saturated  the  first  and  last  arcs  of  an  unsaturated  elementary 
path  are  necessarily  w*  and  w**.  Hence  must  be  such  that  the  flow 
in  arcs  w*  and  w**  is  increased  by  one,  and  arcs  of  types  u*  and  u**  are 
still  saturated.  Thus  is  also  transformable.  If  f.^^  is  not  maximal 
the  argument  may  be  repeated,  so  that  ultimately  ideal  flow  patterns  of 
values!,  i+l,  .  .  .,  M,  which  are  all  transformable,  can  be  generated. 
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From  Lemma  19  it  follows  that  no  feasible  but  nontransformable  flow 
patterns  of  values  i,  i+ 1,  .  .  M  can  be  ideal.  The  proof  is  complete 
if  m  is  taken  as  the  smallest  integer  i  such  that  a  transformable  flow 
pattern  of  value  1  exists. 

Corollary 

Since  there  is  a  feasible  flow  pattern  of  value  k  in  (X,  U)  if 

and  only  if  there  is  a  transformable  x^-»x^  flow  pattern  of  value 
k  +  Sgb(u)  in  (X,  U)^t  follows  that  there  are  feasible  flow  patterns 

of  value  k  for  m  <  k  <  M,  where  m  =  m  -  2gb(u)  and  M  =  M  -  2gb(u), 
and  for  no  other  values. 

PROCEDURE  FOR  SOLVING  PROBLEM  B 

Every  Type  B  problem  can  be  solved  by  the  following  method; 

(a)  Given  a  Type  B  network  (X,  U);  construct  the  associated  Type  A 
network  (X,  U), 

(b)  Solve  the  Type  A  problem  for  this  network,  obtaining  a  sequence 
{  L ,  f„,  .  .  . ,  fr^}  where  each  f.  is  an  ideal  x,  -►x  flow  pattern  of  value  i 
in  (X,  U)  and  M  is  the  Xj^-*x^  capacity  of  (X,  U), 

(c)  If  none  of  the  L’s  are  transformable  there  are  no  feasible  flows  in 
(X,U). 

(d)  If  at  least  one  L  is  transformable  then  for  some  m  it  will  be  true 

that  (7, . f —  ,  }  are  nontransformable  and  {f— - . f-rr}  are 

1  m-1  m  M 

transformable.  It  can  then  be  asserted  that  m  =  m  -  £gb(u)  and 
M  -  M  -  Eg  b(  u )  are  the  smallest  and  largest  values  for  which  feasible 
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X, -•x  flow  patterns  in  (X,  U)  exist.  Moreover  f.  =  (i>f .  .  _  ,  ,  ,  is  an 
1  n  ^  '  1^1+  Sgb(u) 

ideal  flow  pattern  of  value  i  for  i  =  m,  m  + 1,  .  .  . ,  M. 


This  procedure,  incorporating  in  step  b  the  algorithm  of  the  previous 
section,  or  any  other  known  technique  for  solving  Problem  A,  constitutes 


the  complete  solution  to  Problem  B. 
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SOLUTION  OF  PROBLEM  C 


Let  (X,  U)  be  a  weighted,  capacitated  network,  and  let  f  be  a  feasible 
flow  pattern  in  (X,  U).  Suppose  that  one  wishes  to  determine 
whether  or  not  f  is  ideal.  If  the  network  is  a  Type  A  network,  a  theoret¬ 
ical  criterion  has  already  been  given:  f  is  ideal  if  and  only  if  e(C;f )  >  0 
for  every  elementary  cycle  C  in  (X,  W)  where  W  denotes  U  U  U'  as  usual. 
Suppose  f  is  not  ideal.  Since  it  has  been  established  that 

T(f+  g^)  =  T(f)+  e(C;f)  , 

where  g^  is  the  elementary  cycle  flow  associated  with  C,  fj^  =  f  +  g^  is  a 
feasible  flow  pattern  having  the  same  value  as  f  but  reduced  cost,  if 
e(C;f)  <  0.  The  reasoning  can  be  repeated,  using  f^  in  place  of  f.  Either 
e(C;fj^ ) ^0  for  all  C  or  else  one  can  find  C  such  that  ^2  ”  ^ 

smaller  total  cost.  Proceeding  in  this  fashion  an  f^  is  ultimately  obtained 
such  that  e(C;fp)  >  0  for  all  C.  Then  f^  must  be  ideal.  An  algorithm  that 
locates  a  simple  cycle  C  such  that  e{C;f)  <  0  whenever  such  a  cycle  exists 
is  described  in  this  section.  This  algorithm  is  essentially  the  same  as 
Algorithm  1,  which  was  used  earlier  to  solve  Problem  A. 
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Although  Problem  C  is  posed  for  networks  of  Type  B  as  well  as  for 
those  of  Type  A  it  is  sufficient  to  produce  a  general  solution  for  Type  A 
networks.  For  if  (X,  U)  is  a  Type  B  network  let  (X,  U)  denote  the  asso¬ 
ciated  Type  A  network  defined  in  the  previous  section.  Using  the  notation 
of  that  section,  if  f  is  a  feasible  flow  pattern  in  (X,  U)  then  f  =  f 

is  a  transformable  flow  pattern  in  (X,  U),  Moreover  f  is  ideal  in 

(X,  U)  if  and  only  if  f  is  ideal  in  (X,  U).  Hence  f  can  be  tested  to  determine 
whether  or  not  it  is  ideal.  If  it  is  not,  suppose  that  it  is  changed  into  an 
ideal  transformable  pattern  f^^  having  the  same  value.  Then  f^^  =  (t)f^  is  an 
ideal  pattern  in  (X,  U),  having  the  same  value  as  f.  Hence  Problem  C  in 
a  Type  B  network  can  always  be  solved  by  solving  a  related  problem  in  the 
associated  Type  A  network. 

For  Type  A  networks,  the  key  to  solving  Problem  C  is  the  following: 
Theorem  21 

If  Algorithm  1  is  applied  to  a  feasible,  nonideal  flow  pattern  f 

and  if  every  vertex  eventually  has  a  finite  label  then  the  algorithm  will 
continue  to  relabel  vertices  after  n  iterations.  If  y  is  relabeled  during 
the  m^^  iteration,  where  m  >  n,  and  if  the  procedure  is  interrupted 
immediately  after  this  relabeling,  then  the  approach  set  A(y)  will  contain 
a  cycle. 

Proof.  The  second  part  of  the  theorem  will  be  proved  first.  Suppose 
that  y  is  relabeled  on  the  m^**  iteration,  where  m  >  n,  and  the  algorithm  is 
then  halted.  If  A(y)  were  a  simple  Xj^“*y  path  P  at  this  time  it  would  follow 
from  Lemma  14  that  V(y )  >  e( P;f ),  But  P  can  have  at  most  (n  -  1 )  arcs, 
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th 

and  hence  it  follows  from  Lemma  10  that  V(y)  <  e(P;f)  before  the  m^ 

iteration  started,  which  contradicts  the  fact  that  it  was  relabeled  on  the 

iteration.  Hence  A(y)  cannot  be  a  simple  Xj^-»y  path.  It  was  seen 

earlier  that  if  A(y )  was  not  a  path  then  it  contained  a  cycle.  So  the  second 

half  of  the  theorem  has  been  established. 

It  remains  to  show  that  some  vertices  will  in  fact  be  relabeled  on 
th 

iterations  after  the  n*'”.  Since  f  is  not  ideal  there  is  a  simple  cycle  C 
such  that  e(C;f)  <  0.  This  follows  from  Lemma  8.  Suppose  this  cycle  is 

{wj,  Wg . where  Wj^  s  (y^.yg),  ^2~  . 

Wfc  =  Consider  any  time  during  the  algorithm  at  which  all 

vertices  of  the  set  {y^^,  y^,  ....  yj^}  have  finite  labels.  It  has  been 
assumed  that  from  some  point  in  time  onward  this  will  be  the  case. 
Suppose  that  at  such  a  time  the  algorithm  is  halted.  The  following 
inequalities  cannot  all  be  true: 

V(y2)<V(yj)  +e(Wj;f) 

V(y3)<V(y2)  +e(W2;f) 


V(y^)iV(y^.l)  +  e(w^_j;t) 

V(yj);V(y^)  +e(W|j;f) 

for  this  would  imply  that  0  <  e(C;f).  So  at  least  one  label  can  be 
improved.  But  this  will  always  be  the  case,  since  the  algorithm  was 
halted  at  an  arbitrary  time  after  all  vertices  of  the  set  (y^^,  yg,  .  .  . ,  yj^} 
attained  finite  labels.  So  the  algorithm  will  never  arrive  at  a  "  stable" 
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set  of  labels  such  that  V(y)  >  V(x)  +  e(w)  for  every  we  W  where  w  (x,  y). 
This  completes  the  proof. 

Let  f  be  a  feasible  Xj^-*x^  flow  pattern  in  a  Type  A  network  (X,  U), 
with  U'  and  W  defined  as  usual.  A  new  network  (X,  U)  will  now  be  intro¬ 
duced  by  adding  additional  arcs  to  U.  Specifically,  for  every  x^  4=  x^^  an  arc 
V.  is  created,  with  v^  e  (Xj^,x. ).  Flow  bounds  and  costs  are  assigned  as 
follows: 

b(v.)  =0 
c(v.)  =1 
a^(v.)  =  K 

a2(Vi)  =  0. 

Here  K  is  chosen  to  be  an  integer  larger  than  S^max  I  aj^(u),  a2(u)]. 

Thus  K  is  larger  than  the  effective  cost  ( relative  to  f )  of  any  simple  path 
joining  two  vertices  by  means  of  arcs  of  W,  provided  that  none  of  the 
arcs  is  saturated. 

If  Udenotes  U  augmented  by \v„,  v  ,  .  .  .,  v  )andifb(v.),  c(v.), 

^  •t  n  11 

rs/ 

a.  ( V. ),  and  a„(  v. )  are  defined  in  the  above  manner  then  (X,  U)  is  again 

XI  A  1 

a  Type  A  network.  Moreover  if  we  extend  the  definition  of  f  by  setting 
f(u)  =  0  for  u  =  V.  and  i  =  2,  3,  .  .  . ,  n  then  the  extended  function  f  is  a 
feasible  flow  pattern  in  (X,  U).  Ideal  patterns  in  (X,  U)  and  (X,  U) 

are  related  by  the  following  result: 

Lemma  22 

If  f  is  a  feasible  Xj^-»x^  flow  pattern  in  a  Type  A  network  (X,  U)  and  if 
(X,  U)  and  f  are  defined  in  the  indicated  manner  then  f  is. ideal  in  (X, U) 
if  and  only  if  f  is  ideal  in  (X,  U). 
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Proof.  Suppose  that  f  is  a  feasible  flow  pattern  in  (X,  U).  Let  g  be 


another  feasible  flow  pattern  in  (X,  U),  having  the  same  value  as  f,  and 
let  f  and  g  be  these  patterns  extended  to  (X,  U)  by  setting  f(v^)  =  g(v. )  =  0 
for  i  =  2,  3,  .  .  . ,  n.  The  cost  of  f  relative  to  (X,  U)  is  the  same  as  that 
of  f  relative  to  (X,  U)  since  f(u)  =  0  for  ue  (U  -  U).  Similarly  for  g  and 
g.  So  if  f  is  ideal  in  (X,  U)  then  f  is  ideal  in  (X,  U).  For  otherwise  there 
would  be  a  g  in  (X,  U)  with  cost  smaller  than  f,  and  hence  gwith  cost 
smaller  than  f,  which  cannot  be  if  f  is  ideal. 

It  remains  to  show  that  if  f  is  ideal  in  (X,  U)  then  f  is  ideal  in  (X,  U). 
The  argument  here  is  slightly  more  involved,  since  one  must  consider  all 
feasible  flows  of  a  given  value  in  (X,U),  not  merely  those  such  that 


^(Vg)  =  f(V3)  =  .  .  .  =  f(v^)  =  0  , 

in  order  to  determine  an  ideal  flow.  Suppose  that  f  is  ideal  in  (X,  U),  but 
that  f  is  not  ideal  in  (X,  U),  Let  h  be  ideal  in  (X,  U)  and  have  the  same 
value  as  T.  Note  first  that  h  must  be  such  that  h(v. )  *  1  for  some 
i  =  2,  3,  .  .  . ,  n.  For  if  h  is  the  extension  g  of  some  pattern  g  in  (X,  U), 
since  g  and  g  have  the  same  total  cost,  it  may  be  concluded  that  g  is 
"cheaper"  than  f,  since  f  and  f  have  the  same  cost  and  g  is  cheaper  than 
7.  But  this  contradicts  the  assumption  that  f  is  ideal.  Hence  if  f  is  ideal 
but  f  is  not,  then  an  ideal  pattern  h  in  (X,  U)  must  be  such  that  h(  vp  =  1 

rs/ 

for  some  v..  For  convenience  assume  that  h(v„ )  =  1.  Now  f  and  h  are 

1  M 

rsf 

feasible  patterns  of  the  same  value  in  (X,  U).  The  zero-valued  pattern 


f  -  h  can  be  decomposed  into  conformal  simple  cycle  flows  gj^,  gg,  .  .  . , 
One  of  these,  say  g.,  is  such  that  gj(V2)  ”  since  f(V2)  =  0  and  h(V2)  = 


1. 
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From  Lemma  2  it  is  known  that  h  +  g.  is  feasible.  Hence  e(C.;h)  <  «o, 
from  Lemma  6,  where  is  the  simple  cycle  in  W corresponding  to  g^. 

The  cycle  C.  consists  of  some  arcs  of  form  uc  U  and  u'  e  U'  and  of  the 
arc  v^.  For  each  arc  of  form  u  or  u'  the  effective  cost  relative  to  h  is 
clearly  at  most  max  (  aj^(u),  a2(u)}.  For  v^  on  the  other  hand 
e(v^;h)  =  -K  where  K  is  the  unit  cost  aj^(V2).  But  K  >  S^max  {aj^(u), 
a  (u)},  and  hence  clearly  e(C.;h)  <  0.  This  contradicts  the  fact  that  h  is 
ideal.  So  if  f  is  ideal  in  (X,  U),  f  is  ideal  in  (X,  U).  This  completes  the 
proof. 

Note  that  if  f  is  a  feasible  Xj^-*x^  flow  pattern  in  a  Type  A  network 

fsj  /sj  f»a 

(X,  U)  and  if  (X,  U)  and  f  are  defined  in  this  way  then  f  is  feasible  and 
nonmaximal  in  (X,  U).  The  latter  property  follows  from  the  fact  that  we 
could  augment  f  by  a  unit  flow  in  arc  v^^  =  (x^,  x^)  and  thus  increase  the 
value  of  flow  by  one.  Moreover  f  has  the  property  that,  after  one  iteration 
of  Algorithm  1,  every  vertex  of  (X,  U)  will  have  a  finite  label.  This 
follows  from  the  fact  that  e(v.;f)  =  K  <  «»for  i  =  2,  3,  ,  .  . ,  n,  and  that 
directly  connects  Xj^  and  x..  Hence  f  meets  the  requirements  for 
satisfying  either  Theorem  15  or  21,  depending  on  whether  or  not  it  is 
ideal.  These  two  theorems  together  with  Lemma  12  justify  the  following 
general  procedure: 

PROCEDURE  FOR  SOLVING  PROBLEM  C  IN  TYPE  A  NETWORKS 

(a)  Given  a  feasible  Xj^-*x^  flow  pattern  in  (X,  U)  define  (X,  U)  and 
f  in  the  manner  indicated  earlier. 
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(b)  Apply  Algorithm  1  to  f  in  (X,  U)  and  proceed  until  one  of  the 
following  occurs: 

(i)  The  algorithm  terminates,  with  V(y)  <  V(x)  +  e(w;  f)  for 
all  w€  W  where  ws  (x,y). 

(ii)  Vertex  x^^  is  relabeled,  so  that  V(Xj^)  <  0. 

(iii)  (n+  1)  complete  iterations  have  been  accomplished. 

(c)  If  case  (i)  occurs  f  is  ideal;  hence  (Lemma  22)  f  is  ideal.  If 

case  (ii)  occurs  A(Xj^)  contains  a  cycle  C  such  that  e(C;T)  <  0.  Replace 

f  by  the  less  costly  flow  f  +  g  ,  where  g  is  the  elementary  flow  pattern 

c  c 

corresponding  to  C.  If  case  (iii)  occurs,  and  y  is  a  vertex  relabeled 

st  ^ 

during  the  (n+  1)  iteration,  then  A(y)  contains  a  cycle  C  with  e(C;f)  <0. 

Replace  f  by  f.+  g  .  In  either  of  the  last  two  cases  repeat  step  b  with 
c 

^  r*j 

f  +  g  in  place  of  f.  Continue  this  process  until  a  point  is  reached  such 
o 

that  the  outcome  of  step  b  is  b(i).  When  this  occurs  an  ideal  pattern  in 
(X,  U)  has  been  found,  and  also  a  corresponding  ideal  pattern  in  (X,  U). 

The  reason  for  distinguishing  cases  ( ii )  and  ( iii )  in  step  b  of  the 
procedure  is  simply  to  shorten  the  computations  whenever  V(  x^^ )  <  <» 
before  n+ 1  iterations  have  been  completed.  The  procedure  can  be  further 
shortened  by  a  device  essentially  the  same  as  that  employed  in  Lemma  7. 

rs/ 

Once  a  cycle  C  has  been  found  such  that  e(C;  f )  <  0,  evaluate  Q  for  all 

wc  C  as  in  Lemma  7.  Let  Q  =  min  Q  .  Then  f  +  Q  •  g  ,  will  be  a  feasible 

wcC  ^  c 

pattern  whose  cost  is  -Q  •  e(C;f )  less  than  that  of  f.  This  can  be  estab¬ 
lished  rigorously  by  a  line  of  reasoning  very  similar  to  that  employed  in 
the  proof  of  Lemma  7.  ( The  main  point  in  the  reasoning  is  that  Q  is  the 
largest  integer  such  that  one  can  add  up  to  units  of  flow  to  f(w)  before 
the  effective  cost  of  w  changes. ) 
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GLOSSARY  OF  PRINCIPAL  SYMBOLS 


Symbol 

Brief  Explanation 

Page 

X,  y,  z 

Individual  vertices 

7 

X 

Set  of  vertices 

7 

U,  V,  w 

Individual  arcs 

7 

U 

Set  of  arcs 

7 

us(x,y) 

"  u  is  directed  from  x  to  y  ” 

7 

P.  c 

Path  (cycle) 

8 

u' 

"  Reverse"  of  arc  u 

9 

U' 

Set  of  all  u*  for  u  €  U 

9 

W 

UuU* 

9 

f.  g.  h 

Flow  patterns 

11 

f(u),  g(u),  h(u) 

Arc  flows 

11 

U(x-),  Ue*x) 

Set  of  arcs  "  leaving"  f  entering" )  x 

11 

Of(x) 

Net  output  at  x  relative  to  f 

12 

f+gi  f-g 

Arc-by-arc  sum  (difference)  of  f  and  g 

14 

gp.  gQ 

Elementary  path  ( cycle )  flow 

15 

b(u),  c(u) 

Lower  (upper)  bounds  on  flow 

18 

a^(u),  a^iu) 

Unit  cost  in  direction  of  u  ( u’ ) 

20 

Tl  f(u)l 

Cost  of  f  in  u 

20 

T(f) 

Total  cost  of  f:  2  T[  f(u)  ] 

20 

e(u;f) 

Effective  cost  of  arc  relative  to  f 

39 

e(P;f),  e(C;f) 

Effective  cost  of  path  (cycle)  relative  to  f 

40 

V(x) 

"  Label"  of  x 

47 

A(x) 

Set  of  approach  arcs  associated  with  x 

50 

A 

Set  of  all  approach  arcs 

50 

(X.U) 

Augmented  network  for  Problem  B 

58 

(X,U) 

Augmented  network  for  Problem  C 
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INDEX 


Approach 

arc  47 

set  50 

Arc  7 

adjacency  8 

approach  47 

end-points  of  8 

parallelism  8 

strict  adjacency  8 

strict  parallelism  8 

Cost 

effective  39 

of  arc  flow  20 

of  flow  pattern  20 

unit  20 

Cycle  9 

effective  cost  of  40 

flow  15 

saturated  35 

simple  9 

Flow  11 

arc  11 

characteristic  39 

elementary  cycle  15 

elementary  path  15 

feasible  arc  18 

Flow  patterns  11 

conformal  14 

cost  of  20 

difference  of  14 

elementary  15 

feasible  18 

ideal  20 

sum  of  14 

transformable  60 

X— y  14 


value  of 

14 

Graph 

connected 

11 

directed 

7 

finite 

7 

geometric 

8 

Incidence 

8 

Lower  bound  on  flow 

18 

Net  output 

12 

Network 

11 

capacitated 

18 

capacity  of 

24 

type  A 

23 

type  B 

24 

wei^ted 

20 

Path 

8 

effective  cost  of 

40 

flow 

15 

simple 

9 

saturated 

35 

X— y 

8 

Sink 

12 

Source 

12 

Upper  bound  on  flow 

18 

Vertex 

7 

initial 

7 

intermediate 

12 

isolated 

8 

net  ou^t  at 

12 

terminal 

7 
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